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It is shown that a first-order cosmological perturbation theory for Friedmann-LemaJtre- Robertson- 
Walker universes admits one and only one gauge-invariant variable which describes the perturbation 
to the energy density and which becomes equal to the usual energy density of the Newtonian theory 
of gravity in the limit that all particle velocities are negligible with respect to the speed of light. 
The same holds true for the perturbation to the particle number density. 

A cosmological perturbation theory based on these particular gauge-invariant quantities is more 
precise than any earlier first-order perturbation theory. In particular, it explains star formation in 
a satisfactory way, even in the absence of cold dark matter. In a baryon-only universe, the earliest 
stars, the so-called Population III stars, are found to have masses between 400 and 100,000 solar 
masses with a peak around 3400 solar masses. If cold dark matter is present then the star masses 
are between 130 and 13,000 solar masses with a peak around 450 solar masses. They come into 
existence between 100 Myr and 1000 Myr. At much later times, star formation is possible only in 
high density regions, for example within galaxies. Late time stars may have much smaller masses 
than early stars. The smallest stars that can be formed have masses of 0.2-0.8 solar mass, depending 
on the initial internal relative pressure perturbation. 

It is demonstrated that the Newtonian theory of gravity cannot be used to study the evolution of 
cosmological density perturbations. 
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I. INTRODUCTION 



It is well-known, or formulated more precisely, it is generally accepted, that in a universe filled with only 'ordinary 
matter,' i.e. elementary particles and photons but not cold dark matter (cdm), the linear perturbation theory predicts 
a too small growth rate to account for star formation in the universe. In this article we establish that this is not 
true. The reason brought forward in all former treatises on first-order cosmological density perturbations is that the 
growth of a relative density perturbation in the era after decoupling of radiation and matter, given by 



is insufficient for relative density perturbations as small as the observed value (5(idcc) ~ 10^'^ to reach the non-linear 
phase for times t < tp, where tp — 13.7 Gyr, the present age of the universe, and idee = 380 kyr, the time of decoupling 
of matter and radiation. This generally accepted conclusion is suggested by the solutions (E8), with u; = 0, of the 
standard relativistic evolution equation (E6) for linear density perturbations in the era after decoupling of matter 
and radiation: this equation yields a growth rate which is much too low to allow for star formation within 13.7 Gyr. 
Therefore, researchers in the field of structure formation have to assume that a significant amount of cdm had 
contracted already before decoupling in order to explain in their simulations the formation of large-scale structure 
after decoupling [1, 2]. 

The purpose of this article is to show that the formation of structure can be explained whether or not cdm is 
present. Our treatise is independent of a particular system of reference and yields results which describe the evolution 
of small density perturbations in the radiation-dominated era and in the era after decoupling of matter and radiation. 
These remarkable and most satisfactory results are a direct consequence of two facts only. Firstly, we use gauge- 
invariant expressions for the first-order perturbations to the energy density, efl^, and particle number density, nfly 
Secondly, in the dark ages of the universe (i.e. the epoch between decoupling and the ignition of the first stars) a 
density perturbation from which stars will eventually be formed, has to cool down [3, 4] in its linear phase in order to 
grow. Consequently, the growth of density perturbations can only be described by a realistic equation of state for the 
pressure in combination with the combined First and Second Laws of thermodynamics (167). Therefore, we use an 
equation of state for the pressure of the form p — p{n, e), where n is the particle number density and e is the energy 
density. 

Our more sophisticated treatise, which in first-order perturbation theory is explicitly gauge-invariant, and which 
uses an equation of state of the form p ~ p(n,e) rather than of the form p = p(e); makes it possible to explain 
structure formation even in the absence of cdm. 



In this article no assumptions or approximations (other than linearization of the Einstein equations and conservation 
laws) have been made in order to reach our conclusions. The only assumption we have made is that Einstein's General 
Theory of Relativity is the correct theory that describes gravitation in our universe on all scales and from the onset 
of the radiation-dominated era up to the present time. 




(1) 



II. MAIN RESULTS 
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In order to study structure formation in the universe, one needs the linearized Einstein equations. The derivation 
of the evolution equations (201) for relative density fluctuations in a Friedmann-Lemaltre-Robertson- Walker (flrw) 
universe filled with a perfect fluid with an equation of state p — p{n,e), is one of the main subjects of this article. 
Since we use a more general equation of state we are forced to derive all basic equations from scratch, instead of taking 
them from well-known textbooks or renowned articles. This has the advantage that our article is self-contained and 
that our results can easily be checked. 

In Section X A we show, using the background (i.e. zeroth-order) Einstein equations and conservation laws connected 
with the combined First and Second Laws of thermodynamics, that the universe as a whole expands adiabatically. 
This is a well-known result. One of the new results of our treatise is that local density perturbations evolve diabatically. 
This has been made clear in Section XD. Only in the non-relativistic limit, where e = e(n) and p = 0, local density 
perturbations evolve adiabatically. 

In the literature about the subject, all efforts to construct a gauge-invariant cosmological perturbation theory 
yield a second-order differential equation for the density contrast function S = £(i)/e(o), with or without an entropy 
related source term. For that matter, our treatise is no exception, as (201a) demonstrates. In contrast to the 
perturbation theories developed in the literature, we find also a first-order differential equation, namely (201b). This 
equation results from the incorporation of the equation of state p = p{n, e) for the pressure and the particle number 
conservation law {nu'^)-^ — 0, (45). The consequences of equation (201b) are radical: this equation implies that 
density perturbations in the total energy density are gravitationally coupled to density perturbations in the particle 
number density. This is the case for ordinary matter as well as CDM throughout the history of the universe from the 
onset of the radiation-dominated era until the present. As a consequence, perturbations in CDM evolve gravitationally 
in exactly the same way as perturbations in ordinary matter do. The assumption that CDM would have clustered 
already before decoupling and thus would have formed seeds for baryon contraction after decoupling is, therefore, 
questionable. This conclusion has, on different grounds, also been reached by Nieuwenhuizen et al. [5]. This may rule 
out CDM as a means to facilitate structure formation in the universe. 



A. Manifestly Gauge-invariant Perturbation Theory and its Non-Relativistic Limit 

In order to solve the structure formation problem of cosmology, we first develop in Sections IV-XI a manifestly 
gauge-invariant perturbation theory, i.e. both the evolution equations and their solutions are independent of the choice 
of a system of reference. In Section IV we show that there exist two and only two unique and gauge-invariant first-order 
quantities e^, and nfl^ for the perturbations to the energy density and the particle number density. The evolution 
equations for the contrast functions (5^ = efly/e^o) E^nd (5„ = nfl^/n^o^ are given by (201). From their derivation it 
follows that these equations include the — background as well as first-order — Goo- and Goi-constraint equations; 
the Gy-evolution equations; the conservation laws for the energy density; for the particle number density and for 
the momentum; and, finally, the combined First and Second Laws of thermodynamics. Taking into account the 
applicability of equations (201) to the open closed and fiat flrw universes filled with a perfect fiuid described by an 
equation of state p = p{n, e), they are rather simple. 

In Section XII we show that in the non-relativistic limit u/c — > the quantities efl^ and nfl^ survive, and coincide 
with the usual energy density and particle number density [see (229) and (230)], whereas their gauge dependent 
counterparts e^if and n(i) (which are also gauge dependent in the non-relativistic limit) disappear completely from the 
scene. Finally, we show that, in first-order, the global expansion of the universe is not affected by local perturbations 
in the energy and particle number densities. 

B. Large-Scale Perturbations: Confirmation of the Standard Knowledge 

In order to compare our treatise on cosmological density perturbations with the standard knowledge, we consider 
a flat FLRW universe in the radiation-dominated era and in the era after decoupling of matter and radiation. For 
large-scale perturbations these two cases have been thoroughly studied by a large number of researchers from 1946 up 
till now, using the full set of linearized Einstein equations and conservation laws. Consequently, our reflnement of their 
work cannot be expected to give results that differ much from those of the standard theory. Indeed, we have found 
that for large-scale perturbations our manifestly gauge-invariant treatise corroborates the outcomes of the standard 
theory in the large-scale limit of both eras, with the exception that we do not find, of course, the non-physical gauge 
mode <5gaugc oc t~^ which plagues the standard theory. 

For example, our perturbation theory yields in the radiation-dominated era the well-known solutions (258) oc t 
and oc t^^^ [6-11] and in the era after decoupling of matter and radiation we get the well-known solutions (310) 
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5^ cx t^^^ [6-11] and (xt [12, 13]. A new result is, however, that the solutions (258) as well as (310) follow from 
one second-order differential equation, namely (245a) and (293) respectively. 

C. Small-Scale Perturbations and Star Formation: New Results 

The major difference between our treatise and the standard treatise on the subject lays in the evolution of small-scale 
density perturbations. For a radiation-dominated universe, the standard theory yields oscillating density perturbations 
(352) with a decaying amplitude. In contrast, our theory yields oscillating density perturbations with an increasing 
(260) amplitude. This difference is entirely due to the presence of the spurious gauge modes in the solutions of the 
standard equations, as we will explain in detail in Section XVI A. 

After decoupling of matter and radiation at z = 1091, the results of our treatise and the results found in literature 
differ also considerably. Just as is done in previous researches, we take as equations of state for the energy density 
and the pressure e — nm-ac^ + ^nksT and p — nksT, respectively, in the background as well as in the perturbed 
universe. Since mjiC^ 3> fceT throughout the matter-dominated era after decoupling, it follows that one may neglect 
the pressure nk^T and kinetic energy density ^nksT with respect to the rest-mass energy density nmHC^ in the 
unperturbed universe and that in the perturbed universe one has 6^ w (5„. Therefore, one takes in the literature 

= and solves the (homogeneous) second-order evolution equation for S^- As we have shown in Section XIII C 3, 
this yields slowly growing density perturbations with 6p = 6^, i.e. the relative pressure perturbation is equal to the 
relative energy density perturbation, and a vanishing relative matter temperature perturbation St- 

In contrast to the standard method, our perturbation theory yields next to the usual second-order evolution equation 
(293) for also a first-order evolution equation (278) for the difference (5„ — (5^. Therefore, we need not take (5„ exactly 
equal to 5^, so that in our treatise we may have Sp ^ 6^ and St 0. As a consequence, our resulting second-order 
evolution equation (298) becomes inhomogeneous and the initial relative matter temperature perturbation 5T(idcc, x) 
enters the source term. This proves to be crucial for star formation. Although in a linear perturbation theory 
|JT(i,a;)| < 1; this quantity is not constrained to be as small as 6e{tdcc,x) w 5n{tdcc,x) w 10^^, as is demanded 
by WMAP-observations [14-18]. Since the gas pressure p = nk^T is very low, its relative perturbation 5p = pfl)/P(a) 
and, accordingly, the matter temperature perturbation STit^ecx) could be large. We have shown that just after 
decoupling at z = 1091 negative relative matter temperature perturbations as small as —0.5% yields massive stars 
within 13.7 Gyr. The very first stars, the so-called Population ill stars [19-21], come into existence between lO^Myr 
and 10'^ Myr and have masses between 4 x 10^ M0 and 10^ M©, with a peak around 3.4 x 10^ M©. Stars lighter than 
3.4 X 10^ Mq come into existence at later times, because their internal gravity is weaker. On the other hand, stars 
heavier than 3.4 x 10'^ also develop later, since they do not cool down so fast due to their large scale. These 
conclusions, which are valid only in a universe filled with a baryonic fluid, are outlined in Figure 1. However, if CDM 
is present then the peaks in Figure 1 are at 4.5 x 10^ M0, whereas for hot dark matter (hdm) the peaks are found to 
be at 8.2 x 10'^ M©. The relation between the particle mass and the mass of a star will be explained in Section XIV B. 
The peaks in Figure 1 can be considered as the relativistic counterparts of the classical Jeans mass. 

At much later times, star formation is still possible, however the mass of the stars may be much smaller. For example, 
if star formation starts at z = 1 or later then the smallest stars that can be formed have masses of G.2M0-O.8M0, 
depending on the initial internal relative matter temperature perturbation. Also the initial density perturbations 
must be considerable at late times in order to make star formation feasible: for star formation starting at z = 1 
one must have 0.7 < (5„ (5e < 1. This shows that late time star formation is possible only in high density regions 
within galaxies, but not in intergalactic space. These findings are summarized in Figures 2 and 3. In contrast to the 
relativistic theory developed in this article the standard Newtonian theory of linear perturbations, which does not 
follow from Einstein's gravitational theory, predicts, as can be seen from Figure 4, a lower limit for star formation of 
1.7 M0, implying that our Sun could not exist at all! This failure of the standard Newtonian perturbation theory can 
be attributed to the gauge mode which is present in the solution (359) of the standard equation (357). 

We conclude that density perturbations in ordinary matter can account for star formation. There is no need to 
make use of alternative gravitational theories nor of the inclusion of CDM: the Theory of General Relativity can be 
used to explain structure in our universe. The important conclusion must be that Einstein's gravitational theory not 
only describes the global characteristics of the universe, but is also locally successful. 

III. LINEAR PERTURBATIONS IN THE GENERAL THEORY OF RELATIVITY 

Perturbation theories for flrw universes are, in general, constructed along the following lines. First, all quantities 
relevant to the problem are divided into two parts: 'a background part' and 'a perturbation part.' The background 
parts are chosen to satisfy the Einstein equations for an isotropic universe, i.e. one chooses for the background quantities 
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the FLRW-solution. Because of the homogeneity, the background quantities depend on the time coordinate t only. 
The perturbation parts are supposed to be small compared to their background counterparts, and to depend on the 
space-time coordinate x = {ct, x). The background and perturbations are often referred to as 'zeroth-order' and 'first- 
order' quantities respectively and we will use this terminology also in this article. After substituting the sum of the 
zeroth-order and first-order parts of all relevant quantities into the Einstein equations, all terms that are products of 
two or more quantities of first-order are neglected. This procedure leads, by construction, to a set of linear differential 
equations for the quantities of first-order. The solution of this set of linear differential equations is then reduced to a 
standard problem of the theory of ordinary, linear differential equations. 

A. History 

The first systematic study of cosmological density perturbations is due to Lifshitz [22, 23] (1946) and Lifshitz and 
Khalatnikov [6] (1963). They considered small variations in the metric tensor to study density perturbations in the 
radiation-dominated {p — |e) and matter-dominated (p = 0) universe. The use of metric tensor fiuctuations makes 
their method vulnerable to spurious solutions, the so called gauge modes. Adams and Canuto (1975) [7] extended 
the work of Lifshitz to a more general equation of state p = we, where w is a constant. In 1966 Hawking [24] 
presented a perturbation theory which is explicitly coordinate-independent. Instead of using the perturbed metric 
tensor, he considered small variations in the curvature due to density perturbations. In 1976, Olson ]8] corrected 
and further developed the work of Hawking. He defined the density perturbation relative to co-moving proper time. 
The advantage of Olson's method is that the gauge mode, still present in his solutions, can be readily identified 
since gauge modes yield, in his theory, a vanishing curvature perturbation. Bardeen ]12] (1980), was one of the 
first who realized that one should work with variables which are themselves gauge-invariant. He used in his work 
two different definitions of gauge-invariant density perturbations, which, in the small-scale limit, coincide with the 
usual density perturbation which is gauge dependent outside the horizon. Bardeen assumed that a gauge dependent 
perturbation becomes gauge-invariant as soon as the perturbation becomes smaller than the horizon. In Section XII 
we show that this assumption is invalid. Kodama and Sasaki [25[ (1984) elaborated and clarified the pioneering work 
of Bardeen. Ellis, Bruni and Hwang et al. [26, 27[ (1989) and Elhs and van Elst [28[ (1998) criticized the work 
of Hawking and Bardeen and gave an alternative and elegant representation of density fiuctuations. Their method 
is both fully covariant and gauge-invariant. Although the standard equations (355) and (362) are, according to 
Einstein's General Theory of Relativity, closely related via (E6), equation (362) follows from the theory of Ellis et al, 
but equation (355) cannot be derived from their method. Moreover, they did not take into account the perturbed 
constraint equations. Mukhanov, Feldman and Brandenberger, in their 1992 review article [13[ entitled 'Theory of 
Cosmological Perturbations' mentioned or discussed more than 60 articles on the subject, and, thereupon, suggested 
their own approach to the problem. Their method is also discussed in the textbook of Mukhanov [29[. A disadvantage 
of their method is that their perturbation theory does not yield the usual Poisson equation in the non-relativistic 
limit v/c ^ 0, implying that their gauge-invariant quantities cannot be linked to their Newtonian counterparts. For 
an overview of the literature we refer to Mukhanov et al. [13[ and the summer course given by Bertschinger [30[. A 
recent and integral overview of the construction of cosmological perturbation theories for flat FLRW universes is given 
by Malik and Wands [31[. 

The fact that so many studies are devoted to a problem that is nothing but obtaining the solution of a set of 
ordinary, linear differential equations is due to the fact that there are several complicating factors, not regarding the 
mathematics involved, but with respect to the physical interpretation of the solutions. As yet there is no consensus 
about which solution is the best. In this article we will demonstrate that there is one and only one solution to the 
problem how to construct gauge-invariant quantities. This then enables us to solve the problem of structure formation. 

B. Origin of the Interpretation Problem 

At the very moment that one has divided a physical quantity into a zeroth-order and a first-order part, one 
introduces an ambiguity. Let us consider the energy density of the universe, e(a;), and the particle number density of 
the universe, n{x). The linearized Einstein equations contain as known functions the zeroth-order functions s^o)(t) and 
ri(o)(i), which describe the evolution of the background, i.e. they describe the evolution of the unperturbed universe 
and they obey the unperturbed Einstein equations, and as unknown functions the perturbations e^^{x) and n^^{x). 
The latter are the solutions to be obtained from the linearized Einstein equations. The sub-indexes and 1, which 
indicate the order, have been put between round brackets, in order to distinguish them from tensor indices. In all 
calculations, products of a zeroth-order and a first-order quantity are considered to be of first-order, and are retained, 
whereas products of first-order quantities are neglected. 



7 



The ambiguity is that the linearized Einstein equations do not fix the quantities e(i)(a;) and n^i^{x) uniquely. In 
fact, it turns out that next to any solution for ^(i) and of the linearized Einstein equations, there exist solutions 
of the form 

£(i,(x) = + ■0(a;)9oe(o)(t), (2a) 

n(i,(x) = n^r)(x) + ■tp{x)don^o){t), (2b) 

which also satisfy the linearized Einstein equations. Here the symbol do stands for the derivative with respect to 
a;" = ct. The function tpi^) is an arbitrary but 'small' function of the space-time coordinate x = {ct,x), i.e. we 
consider ip{x) to be of first-order. We will derive (2) at a later point in this article; here it is sufficient to note that 
the perturbations e^, and are fixed by the linearized Einstein equations up to terms that are proportional to 
an arbitrary, small function tpi^), usually called a gauge function in this context. Since a physical quantity, i.e. a 
directly measurable property of a system, may not depend on an arbitrary function, the quantities £(1) and n^if cannot 
be interpreted as the real physical values of the perturbations in the energy density or the particle number density. 
But if £(1) and n^^y are not the physical perturbations, what are the real perturbations? This is the notorious 'gauge 
problem' encountered in any treatise on cosmological perturbations. Many different answers to this question can be 
found in the literature, none of which is completely satisfactory; a fact which explains the ongoing discussion on this 
subject. In this article we show that there is a definitive answer to the gauge problem of cosmology. 



IV. GAUGE-INVARIANT QUANTITIES 

In the existing literature on cosmological perturbations, one has attempted to solve the problem that corresponds 
to the gauge dependence of the perturbations £(i) and n(i) (2) in two, essentially different, ways. The first way is 
to impose an extra condition on the gauge field ■^'{x) [30, 32-34]. Another way to get rid of the gauge field ifi{x) 
is to choose linear combinations of the matter variables £(i), and other gauge dependent variables to construct 
gauge-invariant quantities. The latter method is generally considered better than the one where one fixes a gauge, 
because it not only leads to quantities that are independent of an undetermined function, as should be the case for a 
physical quantity, but it also does not rely on any particular choice for the gauge function. 

The newly constructed gauge-invariant quantities are then shown to obey a set of linear equations, not containing 
the gauge function 'il'{x) anymore. These equations follow, by elimination of the gauge dependent quantities in favor 
of the gauge-invariant ones, in a straightforward way from the usual linearized Einstein equations, which did contain 
ip{x). In this way, the theory is no longer plagued by the gauge freedom that is inherent to the original equations 
and their solutions: '4'{x) has disappeared completely, as it should, not with brute force, but as a natural consequence 
of the definitions of the perturbations to the energy and the particle number densities. This method is elaborated 
by Bardeen [12] and Mukhanov et al. [13[. From these two treatises on linear perturbation theory, which differ 
significantly from each other, one is tempted to conclude that gauge-invariant quantities can be constructed in many 
different ways, and that there is no way to tell which of these theories describes the evolution of density perturbations 
correctly. This, however, is not the case, as we will show in this article. 

We follow the method advocated by Mukhanov et al. and Bardeen, but with gauge-invariant quantities which differ 
substantially from those used by these researchers. In fact, we will show that there exist unique gauge-invariant 
quantities 

gi _ ^O^(O) a In \ 

£(1) = £(1) ~ '^TTi ''(I)' I'jaj 
Cof(o) 

"fi) ="(1) - TT^^d)' (3b) 

Cof(o) 

that describe the perturbations to the energy density and particle number density. In these expressions 0(o) and 
are the background and perturbation part of the covariant four-divergence — c^^C/^^^ of the cosmological ffuid 
velocity field U^{x). 

In Section IV A, we will show that the quantities efl^ and nfl^ do not change if we switch from the old coordinates 
to new coordinates x^ according to 

x^^^x'^-^'^ix), (4) 



where the ^^(a;) (/i — 0, 1, 2, 3) are four arbitrary functions, considered to be of first-order, of the old coordinates x^. 
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i.e. (4) is an infinitesimal coordinate transformation, or gauge transformation. In other words, we will show that 

efl^{x) = eflix), (5a) 
hfl{x) = nfl{x), (5b) 

i.e. the perturbations (3) are independent of ^^(x), i.e. gauge-invariant. 

Since the background quantities depend on time, but not on the spatial coordinates, it will turn out that only the 
zero component of the gauge functions ^'^(a:) occurs in the transformation of the first-order gauge dependent variables. 
We will call it iIj{x): 

^{x)^e{x). (6) 

In the perturbation theory, the gauge function 4>{x) is to be treated as a first-order quantity, i.e. as a small (or 
'infinitesimal') change of the coordinates. 

As yet, the quantities (3) are new: they have never been used before. The fact that these quantities are unique 
follows immediately from the linearized Einstein equations for scalar perturbations, and, therefore, cannot be chosen 
arbitrarily. Using the quantities (3), our theory reduces to the usual Newtonian theory (229) and (230) in the limit 
that the spatial part of the cosmological fiuid velocity four-vector [/^ is small compared to the velocity of light. 



A. Construction of Gauge-invariant First-order Perturbations 

We now proceed with the proof that efl^ and nfij are gauge-invariant, i.e. invariant under the general infinitesimal 
coordinate transformation (4). To that end, we start by recalling the defining expression for the Lie derivative of an 
arbitrary tensor field A"'"^ ^...^ with respect to a vector field £,'^{x). It reads 

+ A°' ^ T-'-uV \ti + • • • + A°' ^ ^■■■t£,^ -.y, (7) 

where the semi-colon denotes the covariant derivative. At the right-hand side, there is a term with a plus sign for 
each lower index and a term with a minus sign for each upper index. Recall also, that the covariant derivative in the 
expression for the Lie derivative may be replaced by an ordinary derivative, since the Lie derivative is, by definition, 
independent of the connection. This fact simplifies some of the calculations below. 

Now, let {x''} and {x^^ = x'' — £,^{x)} be two sets of coordinate systems, where £,'^{x) is an arbitrary — but 
infinitesimal, i.e. in this article, of first-order — vector field. Then the components A" ' '^^(...j,(a;) of the tensor A with 
respect to the new coordinates x'' can be related to the components of the tensor A"^'"^ ^^...^{x), defined with respect 
to the old coordinates {x^} with the help of the Lie derivative. Up to and including terms containing first-order 
derivatives one has 

i"-^....(x) = A"-^,...,{x) + {C^Ar-^ ,...,ix) + ■■■. (8) 

For a derivation of this expression, see Weinberg [35], Chapter 10, Section 9. 

Note that x in the left-hand side corresponds to a point, P say, of space-time with coordinates x^ in the coordinate 
frame {x}, while in the right-hand side x corresponds to another point, Q say, with exactly the same coordinates x^, 
but now with respect to the coordinate frame {x}. Thus, (8) is an expression that relates one tensor field A at two 
different points of space-time, points that are related via the relation (4). 

The following observation is crucial. Because of the general covariance of the Einstein equations, they are invariant 
under general coordinate transformations x x and, in particular, under coordinate transformations given by (4). 
Hence, if some tensorial quantity A{x) of rank n {n — 0, 1, . . .) satisfies the Einstein equations with as source term 
the energy-momentum tensor T, the quantity A{x) = A{x) + £^A{x) satisfies the Einstein equations with source 
term T{x) = T{x) + C^T{x), for a universe with precisely the same physical content. Because of the linearity of the 
linearized Einstein equations, a linear combination of any two solutions is also a solution. In particular, C^A, being 
the difference of A and A, is a solution of the linearized Einstein equations with source term C^T. In first-order, 
C^A{x) may be replaced by C^A^o^{t), where A^o){t) is the solution for A{t) of the zeroth-order Einstein equations. 
The freedom to add a term of the form C^A^o){t), with {fi = 0, 1,2,3) four arbitrary functions of first-order, to 
any solution of the Einstein equations of the first-order, is the reason that none of the first-order solutions is uniquely 
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defined, and, hence, does not correspond in a unique way to a measurable property of the universe. This is the 
notorious gauge problem. The additional terms £^A,o){t) are called 'gauge modes.' 
Combining (7) and (8) we have 

A" ^ fj,...^(x) = A"^ ^ fj^...i,{x) A"^ ^ fj.---iy;Ti^ 

+ A— + • • • + A"-\...rC;u. (9) 

We now apply expression (9) to the case that A is a scalar a, a four-vector V^'' and a tensor A^j,^ respectively, 

a{x) = a{x) + C{x)dr(T{x), (10a) 

= + l/M.^^r _ V^^f'.r, (10b) 

For the metric tensor, g^i, we find in particular, from expression (10c), 

9fii^ = 9f,iy + ^ti-.y + iy^i^ (11) 

where we have used that the covariant derivative of the metric vanishes. 

Our construction of gauge-invariant perturbations totally rest upon these expressions for hatted quantities. In case 
a{x) is some scalar quantity obeying the Einstein equations, (t{x) can be divided in the usual way into a zeroth-order 
and a first-order part: 

cr(a;) = (T(o,(t) -t- cr(i)(a;), (12) 

where cr(o)(t) is some background quantity, and hence, not dependent on the spatial coordinates. Then (10a) becomes 

a{x) = a,o,(0 + a^,,{x)+^,°{x)^o(J^,,{t) + e{x)d^a^,,{x). (13) 

The last term, being a product of the first-order quantity S,^{x) and the first-order quantity 9^t(T(i), will be neglected. 
We thus find 

(T(a;) = cr(o)(t) -f (T(i)(a;), (14) 

with 

(T(i)(a;) = cr(i,(a::) 4- ■i/'(a;)9ocr(o)(i), (15) 

where we used (6). Thus, in gauge transformations of scalar quantities, only the zero component of the gauge functions 
need to be taken into account. Similarly, we find from (10b) and (11) 

{x) = vn, + - v^,,e.,r. (16) 

and 

9(i)tLv{x) = g(i)t.u{x) -f + ^u-f,. (17) 

The latter two expressions will be used later. 

We are now in a position that we can conclude the proof of the statement that e^^) and nfl^ are gauge-invariant. 
To that end, we now write down expression (15) once again, for another arbitrary scalar quantity uj{x) obeying the 
Einstein equations. We then find the analogue of expression (15) 

<^(i){x) = uj(i){x) + ilj{x)doUJ^o){t). (18) 

The left-hand sides of (15) and (18) give the value of the perturbation at the point with coordinates x with respect 
to the old coordinate system {a;}; the right-hand sides of (15) and (18) contains quantities with the same values of 
the coordinates, x, but now with respect to the new coordinate system {x}. Eliminating the function ip{x) from 
expressions (15) and (18) yields 

In other words, the particular linear combination occurring in the right-hand side of (19) of any two scalar quantities 
uj and a is gauge-invariant, and, hence, a possible candidate for a physical quantity. 

The expressions (3) are precisely of the form (19). As a consequence, e^^) and n^'j are indeed invariant under the 
general infinitesimal coordinate transformation (4), i.e. they are gauge-invariant. 
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B. Unique Gauge-invariant Density Perturbations 

Expression (19) is the key expression of this article as far as the scalar quantities £(1, and are concerned. It tells 
us how to combine the scalar quantities occurring in the linearized Einstein equations in such a way that they become 
gauge independent. Expression (19) can be used to immediately derive the expressions (3) for the gauge-invariant 
energy and particle number densities. 

In fact, let C/^(x) be the four- velocity of the cosmological fluid. In Section IXC it is shown that in the linear theory 
of cosmological perturbations, defined by the background equations (154) and the perturbation equations (161) only 
three independent scalars play a role, namely 

e{x) = c-^T''''{x)U^{x)U,{x), (20a) 
n{x) ^ c-'^Ni'{x)U^{x), (20b) 
e{x) = c-i[/^;^(x), (20c) 



where 



N^' = nW, (21) 



is the cosmological particle current four- vector normalized according to U^U^ = (? . These scalars are divided 
according to 

£(2:) = e(o)(t) -l-e(i)(a;), (22a) 

n{x) = n(o)(t) + ri(i)(a:), (22b) 

B{x) = B,„,(t) + Q,,,{x), (22c) 

where the background quantities e(o)(t), ?T.(o)(i) and d(o){t) are solutions of the unperturbed Einstein equations. These 
quantities depend on the time coordinate t only. The relation (19) inspires us to consider the gauge-invariant combi- 
nations 

sflix)^eM-^^^u:M, (23a) 

nflix) ^ nM - ^^fM (23b) 

eUx)^6,,,{x)^^^^^u:,,,{x). (23c) 

The question remains what to choose for u in these three cases. In principle, we could choose for uj any of the following 
three scalar functions available in the theory, i.e. we could choose e, n or 0. As we will show in Section XII, the only 
choice which satisfies the perturbed energy density constraint equation (161a) in the non-relativistic limit w/c — >■ is 

uj = e. (24) 

This implies the expressions (3a) and (3b) for the energy and particle number density perturbations, as was to be 
shown. Using (24), we find from (23c) 

<^^a)-§^^a)=0, 0(„^O. (25) 

The physical interpretation of (25) is that, in first-order, the global expansion (20c) is not affected by a local pertur- 
bation in the energy density and particle number density. It should be emphasized here that (25) is not equivalent to 
the 'uniform Hubble constant gauge' of Bardeen [12], i.e. we do not impose the gauge condition 0^^) — 0. In contrast, 
Ofl-^ = follows from the linearized Einstein equations, and, due to its gauge-invariance, holds true in arbitrary systems 
of reference. In other words, a 'uniform Hubble function' is inherent in a relativistic cosmological perturbation theory. 

The fact that the expressions (3) for the gauge-invariant quantities ef^-^ and nfl^ are unique follows immediately 
from the background Einstein equations (154) and their perturbed counterparts (161). Consequently, these quantities 
cannot be chosen arbitrarily. In Section XII on the non-relativistic limit we show that e^i) is the perturbation to the 
energy density and nfi, is the perturbation to the particle number density. 
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V. EINSTEIN EQUATIONS AND CONSERVATION LAWS IN SYNCHRONOUS COORDINATES 

The system of evolution equations (197) or, equivalently (201), the main results of this article, are manifestly 
gauge-invariant. Therefore, one may use any convenient and suitable system of reference to derive these results. 

The choice of a suitable coordinate system can be made as follows. It is well-known that in the Newtonian theory 
of gravity all possible space-time coordinate systems are synchronous, since time and space transformations (220) are 
decoupled in the Newtonian theory. Consequently, in order to show that our perturbation theory yields the Newtonian 
theory of gravity in the non-relativistic limit, it is obligatory to work in a synchronous system of reference. A second 
motivation to use synchronous coordinates is the fact that the background equations (154) are already given with 
respect to synchronous coordinates. Therefore, the evolution equations for scalar perturbations (161) turn out to be 
— in synchronous coordinates — simple extensions of the background equations. This fact has helped us to find the 
scalars (20) which play a key role in the construction of our manifestly gauge-invariant perturbation theory. 

The name synchronous stems from the fact that surfaces with t = constant are surfaces of simultaneity for observers 
at rest with respect to the synchronous coordinates, i.e. observers for which the three space coordinates (z = 
1, 2, 3) remain constant. A synchronous system can be used for an arbitrary space-time manifold, not necessarily a 
homogeneous or homogeneous and isotropic one [23]. 

In a synchronous system of reference the line element for the metric has the form: 

ds^ = c^dt^ - g,j{t, x)dx'dx^. (26) 

In a synchronous system, the coordinate t measures proper time along lines of constant a;'. From (26) we can read 
off that {x° = ct): 

goo{t, x) = 1, go,{t, x) = 0. (27) 

From the form of the line element in four-space (26) it follows that minus gij{t,x), (i = 1,2,3), is the metric of 
three-dimensional subspaces with constant t. Because of (27), knowing the three-geometry in all hyper-surfaces, is 
equivalent to knowing the geometry of space-time. The following abbreviations will prove useful when we rewrite the 
Einstein equations with respect to synchronous coordinates: 

x„- = -i5y, x',=g'''xkj, x'' = +\g'\ (28) 

where a dot denotes differentiation with respect to = d. From (27)-(28) it follows that the connection coefficients 
of (four-dimensional) space-time 

r /^i/ = {9pfj.,i' ~^ 9pv,ti ~ 9piy,p) 1 (29) 
are, in synchronous coordinates, given by 

r%o = Too = r",o = r%, = O, (30a) 
Y'^ij — Kij, F'oj = r*jo = — (30b) 

= Ig''^ {gu,j + gij., - g^J,l) ■ (30c) 

From (30c) it follows that the T'^ij are also the connection coefficients of (three-dimensional) subspaces of constant 
time. 

The Ricci tensor R^j^i, = i?^^Ai^ is, in terms of the connection coefficients, given by 

— r ^u^x r ^A,i/ r ^i/F ^cj f ^^aT ua- (3i) 

Upon substituting (30) into (31) one finds for the components of the Ricci tensor 

Rm^ k'^k- >c\>c'"h (32a) 

Ro^ = - >c''^\k. (32b) 

Rij = - KijK^k + ^Kikti^j + ^Rij, (32c) 

where the vertical bar in (32b) denotes covariant differentiation with respect to the metric of a three-dimensional 
subspace: 

^^j\k = ^^j,k + r*/fc>f'j — r'jfcV;. (33) 
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The quantities ^Rij in (32c) are found to be given by 

3p rk T^fe , T^l j^k T^l j^k tQA\ 

rCij — i ij^k — -L ik,j + i ijJ- kl — i ifci jl- (,o4j 

Hence, ^Rij is the Ricci tensor of the three-dimensional subspaces of constant time. For the components R^\, — g^'^ Rn, 
of the Ricci tensor (32), we get 

R''o^k\-K'k><\, (35a) 

= (35b) 



where we have used expressions (27)-(28). 
The Einstein equations read 



where C^", the Einstein tensor, is given by 



R'j k'j - K\jK''k + ^R'j, (35c) 



G"" = i?^" - iR^c.g^'". (37) 
In (36), A is a positive constant, the well-known cosmological constant. The constant k is given by 

^""^ 2.0766 X 10-43 ni-i kg-^ s^, (38) 



with G = 6.6742 x lO^^m^kg ^ s"^ Newton's gravitational constant and c = 2.99792458 x 10^ ms"^ the speed of 
light. In view of the Bianchi identities one has G^'^-i, = 0, hence, since g^'^-i, — 0, the source term T^'^ of the Einstein 
equations must fulfill the properties 

T^''',^ = 0. (39) 

These equations are the energy-momentum conservation laws. 

In order to derive simultaneously the background and first-order equations, we rewrite the Einstein equations (36) 
in an alternative form, using mixed upper and lower indices: 

R^^^ = k(T^, - i5^,T"„) - A5^,. (40) 

Upon substituting the components (35) into the Einstein equations (40), and eliminating the time derivative of x'^k 
from the ii^o-equation with the help of the i?'j-equations, the Einstein equations can be cast in the form 

- 3i? - x^x'k - 2{kT% + A), (41a) 

^\\^ - = (41b) 

ic'j - x',x\ + ^R'o = KiT'j - \S',T%) - A(5''„ (41c) 

where 

3i? = g'-'3i?,^ =3^^, (42) 

is the curvature scalar of the three-dimensional subspaces of constant time. The (differential) equations (41c) are the 
so-called dynamical Einstein equations: they define the evolution (of the time derivative) of the (spatial part of the) 
metric. The (algebraic) equations (41a) and (41b) are constraint equations: they relate the initial conditions, and, 
once these are satisfied at one time, they are satisfied automatically at all times. 

The right-hand side of equations (41) contain the components of the energy momentum tensor T^^, which, for a 
perfect fiuid, are given by 

T^', = {e+p)u^u,-p5^',, (43) 

where u'^it^x) = c'^U^{t^x) is the hydrodynamic fiuid four-velocity normalized to unity {u^u^ — 1), e{t^x) the 
energy density and p[t, x) the pressure at a point [t, x) in space-time. In this expression we neglect terms containing 
the shear and volume viscosity, and other terms related to irreversible processes. The equation of state for the pressure 

p = p{n,e), (44) 
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where n{t,x) is the particle number density at a point {t,x) in space-time, is supposed to be a given function of n 
and e (see also Appendix A for equations of state in alternative forms) . 

As stated above already, the Einstein equations (41a) and (41b) are constraint equations to the Einstein equations 
(41c) only: they tell us what relations should exist between the initial values of the various unknown functions, in 
order that the Einstein equations be solvable. In the following, we shall suppose that these conditions are satisfied. 
Thus we are left with the nine equations (41c), of which, because of the symmetry of g^j, only six are independent. 
These six equations, together with the four equations (39) constitute a set often equations for the eleven (6 + 3 + 1-1-1) 
independent quantities gij, u*, e and n. The eleventh equation needed to close the system of equations is the particle 
number conservation law N'^^^^ — 0. Using (21), we get 

(nu^);^ = 0, (45) 

where a semicolon denotes covariant differentiation with respect to the metric tensor g^^y. This equation can be 
rewritten in terms of the fluid expansion scalar defined by expression (20c). Using (30), we can rewrite the four- 
divergence (20c) in the form 

6* = u° - x^ku^ + (46) 

where the three-divergence is given by 

^? = w'|fe. (47) 

Using now expressions (20c), (30), (33) and (46), the four energy-momentum conservation laws (39) and the particle 
number conservation law (45) can be rewritten as 

fm ^ j.Qk^^ x^iT\ - >f\.T™ ^ 0, (48a) 

+ ~ 2x\T''" - >f\.r*" = 0, (48b) 

and 

nu° + n^ku'' + nd = 0, (49) 

respectively. Since T^* is a vector and T*-' is a tensor with respect to coordinate transformations in a subspace of 
constant time, and, hence, are tensorial quantities in this three-dimensional subspace, we could use in (48) a bar to 
denote covariant differentiation with respect to the metric gij(t,x) of such a subspace of constant time t. 

The Einstein equations (41) and conservation laws (48) and (49) describe a universe filled with a perfect fluid and 
with a positive cosmological constant. The fluid pressure p is described by an equation of state of the form (44): in 
this stage we only need that it is some function of the particle number density n and the energy density e. 

We have now rewritten the Einstein equations and conservation laws in such a way that one can easily derive the 
background and perturbation equations. 



VI. ZEROTH- AND FIRST-ORDER EQUATIONS FOR THE FLRW UNIVERSE 

We will now limit the discussion to a particular class of universes, namely the collection of universes that, apart from 
a small, local perturbation in space-time, are homogeneous and isotropic, the Friedmann-Lemaitre-Robertson- Walker 
(flrw) universes. 

We expand all quantities Q in the form of series, and derive, recursively, equations for the successive terms of these 
series. We will distinguish the successive terms of a series by a sub-index between brackets: 

Q = (3(0) + vQ(i) + v^Q(2-, + • • • , (50) 

where the sub-index zero refers to quantities of the unperturbed, homogeneous and isotropic FLRW universe. 

In expression (50) i] (77 = 1) is a bookkeeping parameter, the function of which is to enable us in actual calculations 
to easily distinguish between the terms of different orders. 
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A. Zeroth-order Quantities 

This section is concerned with the background or zeroth-order quantities occurring in the Einstein equations. All 
results of this section are standard [35], and given here only to fix the notation unambiguously. 
For a FLRW universe, the background metric g(o)p,iy is given by 



g(omit,x) = 1, g^o)Oiit,x) ^ 0, (51a) 

1 

a^ty 



g(oyij{t,x) = -a'^{t)gij{x), g'(ly^{t,x) = --^^5*^(a;). (51b) 



where gij{x) is the metric of a three-dimensional maximally symmetric subspace: 

5ij = diag ^Y"- ^^,r2,r^sin^n7^ , fc = 0,±l. (52) 

The minus sign in (51b) has been introduced in order to switch from the conventional four-dimensional space-time 
with signature (+, — , — , — ) to the conventional three-dimensional spatial metric with signature (-I-, +, +). 

All background scalars depend on time only. Furthermore, four-vectors have vanishing spatial components. Since 
is a unit vector we have 

<o, - 5%. (53) 

The time derivative of the three-part of the metric 5(o)ij, >f(o)ij, may be expressed in the Hubble function 'H{t) = 
{da / dt) / a{t) . We prefer to use a function 

H(t).«, (54) 
which we will call Hubble function also. Recalling that a dot denotes differentiation with respect to ct, we have 

H=^. (55) 
a 

Substituting the expansion (50) into the definitions (28), we obtain 

X(oyi] = -Hg^o^.j, xl^^j = ~H5''j, = -Hgll^, (56) 

where we considered only terms up to the zeroth-order in the bookkeeping parameter rj. 

Similarly, with (46), (47), (50), (53) and (56) we find for the background fluid expansion scalar, O^^^, and the 
three-divergence, i9(o), 

0,0, = 3i/, = 0. (57) 

These quantities and their first-order counterparts will play and important role in our perturbation theory. 
Using (43), (50), (51) and (53) we find for the components of the energy momentum tensor 

^(0)0 — S(f))! "^(0)0 — 0' ^(o)i = ~P{o)5 ji (58) 

where the background pressure P(o) is given by the equation of state (44), which, for the background pressure, is 
defined by 

P(o) =P(o)("(o),e(o))- (59) 
The background three-dimensional Ricci tensor, (34), is given by 

^R{0)ij = ^{o)ij,k ~ ^(0)ik,j + -'^(0)^^(0)*:/ ^ ^ (o)ik^ (otjl^ (^^) 

where the connection coefficients T^g^^^j are given by 

rfo)ii = ^9(0) {9{0)li,j + 5(0)ij,i — g{0)i],l) , (61) 
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where g\l-^ and 5(o)ij depend on time. Hence, the connection coefficients Fj^^^j^ are equal to the connection coefficients 
y^ij of the metric gij: 

rfo,., = r'=,, ^ \f' {gu,, + gy.. - ~9^,,i) ■ (62) 



Therefore, they do not depend on time. 

Substituting (51) and (52) into (60), combined with (62), we find 



From (63) we have 



%ow^^kg,j. (63) 



2k 

'^,ou{t)^-^f,. (64) 
implying that the zeroth-order curvature scalar ^i?(o) = .9(0) ^R(a)ij is given by 

%0^{t)^-^y (65) 

This quantity and its perturbed counterpart will play an important role in our perturbation theory. Note, that in 
view of our choice of the metric (+, — , — , — ), spaces of positive curvature k have a negative curvature scalar '^i?(o). 
Thus, we have found all background quantities. 



B. First-order Quantities 

In this section we express all quantities occurring in the Einstein equations in terms of zeroth- and first-order 
quantities. The equations of state for the energy and pressure, s{n,T) and p{n,T), are not specified yet. 

Upon substituting the series (50) into the normalization condition w'^w^ = 1, one finds, equating equal powers of 
the bookkeeping parameter rj, 

<) - 0, (66) 
for the first-order perturbation to the four-velocity. Writing the inverse of gki = g{o)ki + Vdwki -t- ■ • • as 

/'=5S+wf/, + --- , (67) 
where g'^j^ is the inverse, oi gio)ki, (51b), we find 



and 



It is convenient to introduce 



so that 



ki _ _ ki U 



9m^ = -9t)9mii- (69) 



Kj = -9(i}ij, (70) 



ff'^g^,, h^,=g^,hu,. (71) 
For the time derivative of the first-order perturbations to the metric, X(i)ij (28), we get 



{1)1] — 2^ij^ -^(l)j ~ 2 i' "^(1) ~ 2^"'- C''^) 



The first-order perturbation 6'(i) to the fluid expansion scalar 6 (46), can be found in the same way. Using (50) and 
(53) one arrives at 

0(1, = - (73) 
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where we used (66) and (72). This expression will play an important role in the derivation of the first-order pertur- 
bation equations in Section VIII. The first-order perturbation iJji) to the three-divergence ?? (47), is 

^iD^^Uk^ (74) 

where we have used that 

(wV)(i) (75) 

which is a consequence of ^'^i^iku\o^ = as follows from (53). 

Upon substituting the series expansion (50), and into (43) and equating equal powers of r/, one finds for the 
first-order perturbation to the energy-momentum tensor 

^a)0 = (76a) 
ra,o = (£(o)+P(o,)<,, (76b) 
Tl,, = -p,„<5'„ (76c) 

where we have used (53) and (66). The first-order perturbation to the pressure is related to e^, and n^i^ by the 
first-order perturbation to the equation of state (44). We have 

P(l) =Pn?l(l) +^££(1), (77) 

where p„ and are the partial derivatives of p{n, e) with respect to n and e. 

Since we consider only first-order quantities, the partial derivatives are functions of the background quantities only, 
i.e. 

Pn = Pn {n^o) , £(0) ) , Pe =Pe (?l(o) , £(o) ) ■ (79) 

Using (30c) and the series (50), we find for the first-order perturbations of the connection coefficients 

rfi)ij — ~9*'io)9(i)imX"o)ij + i^fo) i9(i)ii,j + 9{i)i],i — • (80) 

The first-order perturbations Y^^^ij (80), occurring in the non-tensor ij, happen to be expressible as a tensor. Indeed, 
using (70), one can rewrite (80) in the form 

rfi)y = -i9fo)(^/ib- + - hij\i). (81) 
Using the expansion (50) for ^R^j and V^ij, one finds for the first-order perturbation to the Ricci tensor (34) 

-^(^)ij (l)ij,k ^ {i)ik,j ' (Ojij-*- {l)kl ' ^ {0)kl ^ (ojifc-*- (l)ife-'- (O)ji' \°^) 

which can be rewritten in the compact form 

'^R{i)ii — ^\i)ij\k ~ ^\i)ik\j- (83) 

By substituting (81) into (83), one can express the first-order perturbation to the Ricci tensor of the three-dimensional 
subspace in terms of the perturbation to the metric and its covariant derivatives: 



'^R{i)ij — ~\9^o){hii\'j\k + hiju\k — hijmk — ^/fciih)- (84) 



The perturbation ^R\i)j is given by 

'^^(i)i = (f*^ ^Rpi)m = 5(0) ^R(i)pj + 5 ^Rmh"^], (85) 
where we have used (51), (63), (65) and (71). Upon substituting (84) into (85) we get 

~\9\a){h^ p\j\k + j\p\k - h^k\p\j) + h9mf^^j\k\i + 5 ^R[o)h^j- (86) 
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Taking i = j in (86) and summing over the repeated index, we find for the first-order perturbation to the curvature 
scalar of the three-dimensional spaces 

^^(1) — 9li){h^k\i\j — h^i\]\k) + I ^R{0)h^k- (87) 

This expression will play an important role in the evolution of density perturbations, equations (195), and in the 
non-relativistic limit in Section XII. 

We thus have expressed all quantities occurring in the relevant dynamical equations, i.e. the system of equations 
formed by the Einstein equations combined with the conservation laws, in terms of zeroth- and first-order quantities 
to be solved from these equations. In the Sections VIC and VI D below we derive the background and first-order 
evolution equations respectively. To that end we substitute the series (50) into the Einstein equations (41) and 
conservation laws (48) and (49). By equating the powers of rf , 7]^, . . . , we obtain the zeroth-order, the first-order 
and higher order dynamical equations, constraint equations and conservation laws. We will carry out this scheme for 
the zeroth- and first-order equations only. 



C. Zeroth-order Equations 



With the help of Section VIA and the series (50) we now can find from the Einstein equations (41) and conservation 
laws (48) and (49) the zeroth-order Einstein equations and the conservation laws. Furthermore, in view of the 
symmetry induced by the isotropy, it is possible to switch from the six quantities gij and the six quantities >Cij to the 
curvature '^i?(o)(t) and the Hubble function H{t) only. 



1. Einstein Equations 



Upon substituting (56) and (58) into the (0, 0)-component of the Einstein equations, (41a), one finds 

3ij2 _ i = Ke,o, + A. (88) 

The (0, j)-components of the Einstein equations, (41b), are identically fulfilled, as follows from (56) and (58). We 
thus are left with the six («, j)-components of the Einstein equations, (41c). In view of (56), (58) and (64) we find 

that xlo-fi — >ffo)2 — ^fo)3; "^(0)1 — ^(0)2 = ^(o)3 ^^'^ '^^(^o)i ~ ^^'{0)2 ~ '^^fo)3' whereas for i ^ j these quantities vanish. 
Hence, the six (i, j)-components reduce to one equation, 

H = -3H^ + I + ^^(£(0, ~ P(o,) + A. (89) 

In equations (88) and (89) the background curvature ^R^o) is given by (65). It is, however, of convenience to determine 
this quantity from a differential equation. Eliminating a(t) from expressions (55) and (65) we obtain 

%o,+2H%„,^0, (90) 

where the initial value ^^R^a){tQ) is given by 

fik 

'RM^-^V^y (91) 

in accordance with (65). It should be emphasized that equation (90) is not an Einstein equation, since it is equivalent 
to expression (65). It will be used here as an ancillary relation. 



2. Conservation Laws 



Upon substituting (56) and (58) into the 0-component of the conservation law, (48a), one finds 

£(0, -f3if(e,o, -fp(o,) = 0, (92) 

which is the relativistic background continuity equation. The background momentum conservation laws (i.e. the 
background relativistic Euler equations) are identically satisfied, as follows by substituting (56) and (58) into the 
spatial components of the conservation laws, (48b). 
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The background particle number density conservation law can be found by substituting (53) and (57) into equation 
(49). One gets 

n(o) + 3i?n(„) = 0. (93) 

It can be shown by differentiating the constraint equation (88) with respect to time that the general solution of the 
system (88) and (90)-(93) is also a solution of the dynamical equation (89). Consequently, equation (89) need not be 
considered anymore. This concludes the derivation of the background equations. 



D. First-order Equations 

In this section we derive the first-order perturbation equations from the Einstein equations (41) and conservation 
laws (48) and (49). The procedure is, by now, completely standard. We use the series expansion (50) in rj for the 
various quantities occurring in the Einstein equations and conservation laws of energy-momentum, and we equate the 
coefhcients linear in t] to obtain the 'linearized' or first-order equations. 



1. Einstem Equations 



Using the series expansions (50) for "^i?, k^j and T^o, in the (0, 0)-component of the constraint equation (41a), one 
finds 

2^fo,fc>^(i,/ - '^(1) - 2>^fo)/^\,fc = ^^T%. (94) 

With the zeroth-order expressions (56), the abbreviations (72) and the expression for T^^g, (76), we may rewrite this 
equation in the form 

Hh''k + l%i,=~KS^,y (95) 



Using the series expansion (50) for x'^j and T"i, we find for the (0, i)-components of the constraint equations (41b) 

~ ^a)i\k = (96) 

where we noted that 

(^*i|fc)(i) = (9'^) 

which is a consequence of ^lijik^lojj ~ ^\i)jk^lo)i — which, in turn, is a direct consequence of (56). From (72) 
and (76) we find 

h''k\t - h'',\k = 2K(e(o) + P(o))u(in- (98) 

Finally, we consider the (i, j)-components of the Einstein equations (41c). Using the series expansion (50) for >f*j, 
T% and ^ij^j, we find 

^mi ~ '^(i)j'^(o)fe ~ -^(o)j-^a)fc + '^^(i)j ~ ~ 2^ j-^u)/^)- (99) 

With (56), (72) and (76), we get 

h'j + iHh'j + S'jHh^u + = -K<5^(e(i) -P(i)), (100) 

where ^Rl^j is given by expression (86). 

Note that the first-order equations (95) and (100) are independent of the cosmological constant A: the effect of the 
non-zero cosmological constant is accounted for by the zeroth-order quantities [cf. equations (88) and (89)]. 
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2. Conservation Laws 

We now consider the energy conservation law (48a). With the help of the series expansion (50) for >f*j and T^^, 
one finds for the first-order equation 

^00 I rpOk I k rpl , k rpl k rynOO ,k T^OO n {^C\^\ 

^ (1) + (1) |fc + ^(0)1-^ (i)fc + ^{1)1-^ (0)fc ~ ^(o)k-'-(i) ~ ^(i)k-^ (0) — ^' \^^^) 

where we have used that for a three-vector T'^^ we have 

(T"V)(i,=r°t|fc, (102) 

see expression (75). Employing (56), (58), (72)-(74) and (76) we arrive at the first-order energy conservation law 

£(1) -l-3iJ(e(i) +P(i)) + (£(0) +P(o))0(i) = 0. (103) 

Next, we consider the momentum conservation laws (48b). With the series expansion (50) for j and T^^'^ , we find 
for the first-order momentum conservation law 

T'" + {T'*'\k)m - 2>f(o)fcT'(^i" - 2>f(j)j.T,o" - >C(o)kT'" - x^^k^lo) = 0- (104) 

Using that 

(r'V)a, = ^9%Pi^■,\k, (105) 
and expressions (56), (58), (72) and (76) we arrive at 

\tt ^^^"^ +^'<«))"(i) ~ 9lo)Pw\k + 5-ff(£(o) +P(o))U(i) = 0, (106) 

where we have also used that the covariant derivative of (7*0) vanishes: <7'|Jj|j. — 0. 

Finally, we consider the particle number density conservation law (49). With the series expansion (50) for n, 9, and 
M^', it follows that the first-order equation reads 

?T-(0)U°i) + "-(i)M°o) + »^(0),fcwfi) + J^cij.feufo) + '^(O)^'(i) + "(i)^'(o) = 0. (107) 

With the help of (53), (57) and (66) we find for the first-order particle number conservation law 

n(i, + 3Hn^r) + n^„^0^,^ = 0. (108) 
This concludes the derivation of the basic perturbation equations. 

3. Summary 

In the preceding two Sections VI D 1 and VI D 2 we have found the equations which, basically, describe the pertur- 
bations in a FLRW universe, in first-order approximation. They are equations (95), (98), (100), (103), (106) and (108). 
For convenience we repeat them here 

Constraints: Hh'^k + 5 ^^(i) = —Ke^^, (109a) 

h''k\i - h''i\k = 2k(£(o) -|-p(o))M(i)i, (109b) 

Evolution: h'j + iHh\j + 5\jHh'^k + 2 ^i^j^^. -nS'^ie^^^ - p^,), (109c) 

Conservation: £(1, + 3i7(£(i) +P(i)) + (£(o) +P(o))^(i) = 0, (109d) 
1 d 



c 



dt 



(£(0) +P(o))M(i)J - g\o)Pm\k + 5-ff(£(o) -Hp(o))U(i) = 0, (109e) 
4- 3iln(i)+ 71(0,61(1, =0, (109f) 

where 6'(i,, '^R\i-^j and are given by (73), (86) and (87) respectively. Hence, the equations (109) essentially are 
fifteen equations for the eleven (6 -I- 3 -I- 1 -I- 1) unknown functions h'' j, u^i,, £(i, and n(i,. The pressure P(o, is given by 
an equation of state (59), and the perturbation to the pressure, P(i,, is given by (77). The system of equations is not 
over-determined, however, since the four equations (109a) and (109b) are only conditions on the initial values. These 
initial value conditions are fulfilled for all times t automatically if they are satisfied at some (initial) time t = t^. 
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VII. CLASSIFICATION OF THE SOLUTIONS OF FIRST-ORDER 



It is well-known that, for flat FLRW universes, the set of linear equations (109) can be divided into three independent 
sets of equations, which, together, are equivalent to the original set. In this section we show that this can also be done 
for the open and closed FLRW universes. We will refer to these sets by their usual names of scalar, vector and tensor 
perturbation equations. We will show that the vector and tensor perturbations do not, in first-order, contribute to 
the physical perturbations efl^ and nfly As a consequence, we only need, for our problem, the set of equations which 
are related to the scalar perturbations. By considering only the scalar part of the full set of perturbation equations 
we are able to cast the perturbation equations into a set which is directly related to the physical perturbations efl^ 
and nfly This is the subject of Section VIII. 

At the basis of the replacement of one set (109) by three sets of equations stands a theorem proved by York [36] 
and Stewart [37], which states that a symmetric second rank tensor can be divided into three irreducible pieces, and 
that a vector can be divided into two irreducible pieces. Here, we will use this general theorem to obtain equations 
for the scalar irreducible parts of the tensors h'^j and ^R%)j ^^id the vector u.(i), namely '^R^-^^^^j and 

For the perturbation to the metric, a symmetric second rank tensor we have in particular 

h'^^h^^ + h^^ + K^, (110) 

where, according to the theorem of Stewart ]37], the irreducible constituents /ly^ , h\j and hlj have the properties 



hik^O, (111b) 
*k — 



Kk^O, h'»ik=0, (111c) 



with (j){t,x) and ({t,x) arbitrary functions. The contravariant derivative A^"^ is defined as g^^g-^Ay. The functions h^^^, 

and h^^ correspond to scalar, vector and tensor perturbations respectively. 
In the same way, the perturbation to the Ricci tensor can be decomposed into irreducible components, i.e. 

(112) 



(i)lli (i)-Li ^ 



The tensors "^i?* ,, , ■^i?* and '^^m*, have the properties comparable to (111), i.e. 



'^a,lb = |(^'^^- + '^'''b)' (113a) 
= 0, (113b) 

where ^{t, x) and 7r(i, a;) are two arbitrary functions. By now using (86) for each of the irreducible parts we find 

1 r 



+ C'lfe'^, + C'^l/lfe + i X,('5^-0 + C'^b)] , (114a) 

^^(i)-Li ^ ^ iff(o)(^ip|i|fe + ^^^j\p\k) + \s*'J,)h\j\k\i + I ^R(a)h'^j, (114b) 

= ^ l5(o)(^*p|j|fc + ^%\p\k) + l.9fo)^lj|fe|i + \ ^R{o)h\y (114c) 

Combining expressions (113b) and (114b) it follows that has the property 

h%\i = 0, (115) 

in addition to the property (111b). In Section VII B we show that this additional condition is needed to allow for the 
decomposition (124). 

Combining expressions (113c) and (114c), we find that h\j must obey 

9^^iKX\i\m\l + ^Ti\k\m\l ~ Ki\k\l\m) — 0' (116) 
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in addition to (111c). The relations (116) are, however, fulfilled identically for FLRW universes. This can eas- 
ily be shown. First, we recall the well-known relation that the difference of the covariant derivatives A^'"^ /....^p^q 
and A^'"^ k---i\q\p of an arbitrary tensor can be expressed in terms of the curvature and the tensor itself (Weinberg [35], 
Chapter 6, Section 5) 

^ k---l\p\q ~ ^ k---l\q\p — 

4- , ^ . . . _l , 3ns 



(O)kpq ' ^ k---s ^^{o)lpq 

i 

(O)spq 



k---l^R\o)spq ■■■ A ^ k---l^Rla)spqi (H''') 



where '^R\„-^j^.i is the Riemann tensor for the spaces of constant time. At the right-hand side, there is a term with a 
plus sign for each lower index and a term with a minus sign for each upper index. 
We apply this identity taking for A the second rank tensor h\j to obtain 

l,m _ ini t m 3ps _ us 3nm (^ i q\ 

"'*k\i\m "'*k\m\i ~ "'*s ^'■(a)kim "'*k ^^{o)sim- 

Now note that /i^|„j|j vanishes in view of (111c). Next, we take the covariant derivative of (118) with respect to x', 
and contract with 5'^' 

9'''' ^Tk\i\m\l = 9^\h'^s ^Rto)kim ~ Kk ^R^)sim)\l- (119) 

Next, using the expression which one has for the Riemann tensor of a maximally symmetric three-space, 

'Rto,bcd^k{S\gM~S%~gtc), (120) 
(where k — 0, ±1 is the curvature constant) we find 

5''/irfc|,|™|/ = 0, (121) 
i.e. the first term of (116) vanishes. The second and third term can similarly be expressed in the curvature 

_ _ I m 3ps , um 3ns _ /,s Spm (^00\ 

"'*i\k\m\l "'*i\k\l\m ~ "'*s\k -'^(o)imi "•*i|s ^'■(O)kml "*i|fc ^''(o)smi' K'-'^'^l 

where the general property (117) has been used. Upon substituting the Riemann tensor (120) and contracting with g'^^ 
we then arrive at 

9''\h'^i\k\m\l ~ Kl\k\l\m.) = 0> (123) 

i.e. the second and third term of (116) together vanish. Hence, for FLRW universes, equation (116) is identically 
fulfilled. Consequently, the decomposition (111c) imposes no additional condition on the irreducible part hl^ of the 
perturbation h'j. 

The three-vector can be uniquely divided according to [37] 

■"(1) = "(1)11 +'"(i)-L, (124) 
where it(i)|| is the longitudinal part of with the properties 

V A (tt(i,||) = 0, V-M,i, = V-tt(i,||, (125) 
and M(i)_L is the transverse part of with the properties 

V-M(i,_L=0, VAW(i) = VAW(i,_L, (126) 
where the divergence of the vector tt(i) is defined by [see (47)] 

V-M(i) (127) 
and the rotation of the vector u^J^^ is defined by 

(V A M(i)), EE e/'''u(i)j|fe = e/'''u(i)j- fe, (128) 

where ei^'' is the Levi-Civita tensor with ei^"^ = +1. In expression (128) we could replace the covariant derivative by 
the ordinary partial derivative because of the symmetry of T^jk- 

Having decomposed the tensors h^j, '^R\i-,j and u\^^ in a scalar ||, a vector _L and a tensor part *, we can now 
decompose the set of equations (109) into three independent sets. The recipe is simple: all we have to do is to 
append a sub-index || , _L or * to the relevant tensorial quantities in equations (109). This will be the subject of the 
Sections VII A. VII B and VII C below. 
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A. Tensor Perturbations 



We will show that tensor perturbations are not coupled to, i.e. do not give rise to, density perturbations. Upon 
substituting = hlj and ^Rlr)j — i^ito the perturbation equations (109) and using the properties (111c) 

and (113c), we find from equations (109a), (109b) and (109d) 

£(1) = 0, = 0, 71(1) = 0, = 0, (129) 

where we have also used (77). With (129), equations (109e) and (109f) are identically satisfied. The only surviving 
equation is (109c), which now reads 

K^+3Hhlj+2^Rl^^j=0, (130) 
where ^-Rj^^^ is given by (114c). Using (73), (74), (111c) and (129) it follows from (3) that 

ef;, =0, nfl=0, (131) 

so that tensor perturbations do not, in first-order, contribute to physical energy density and particle number density 
perturbations. Hence, the equations (130) do not play a role in this context, where we are interested in energy density 
and particle number density perturbations only. 

The equations (130) have a wave equation like form with an extra term. Therefore, these tensor perturbations 
are sometimes called gravitational waves. The extra term 3Hhlj in these equations is due to the expansion of the 
universe. The six components /i*^ satisfy the four conditions (111c), leaving us with two independent functions k^j. 
They are related to linearly and circularly polarized waves. 



B. Vector Perturbations 



We will show that, just like tensor perturbations, vector perturbations are not coupled to density perturbations. 
Upon replacing h^j by h\j and ^-R(i)j by ■^i?*ijj^j in the perturbation equations (109), and using the expressions (111b) 
and (113b), we find from equation (109a) and the trace of equation (109c) 

£(1) = 0, P(i, = 0, 77.(1) = 0, (132) 

where we have also used (77). 

Since h^j_j is traceless and raising the index with g^^^ in equation (109b) we get 

/75, + 2^/7^, - 2K(e,o, + P(o,)77fi,, (133) 

where we have used (28) and (56). We now calculate the covariant derivative of equations (133) with respect to , 
and use (115) to obtain 

V • M(i, = 0, (134) 

where we made use of the fact that the time derivative and the covariant derivative commute. With (124)-(126) we 
see that only the transverse part of M(i), namely M(i)_l, plays a role in vector perturbations. From (111b) and (132) it 
follows that the equations (109d) and (109f) are identically satisfied. The only surviving equations are (109b), (109c) 
and (109e), which now read 

= -2K(e(o) + P(o))m(i)_Lj, (135a) 
/il^ + SHh^j + 2 = 0, (135b) 

+ 5iJ(£(o, + P(o))wji,± = 0, (135c) 

where ^Rlij±j is given by (114b). 

Using (73), (74), (132) and (134) we get from (3) 

sfl=0, nfl=0, (136) 



i a 
cdf. 



(£(0) +P(0))W(i)J 
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implying that also vector perturbations do not, in first-order, contribute to physical energy density and particle 
number density perturbations. Hence, the equations (135) do not play a role when we are interested in energy density 
and particle number density perturbations, as we are here. Vector perturbations are also called vortices. 

Since vector perturbations obey V • — 0, they have two degrees of freedom. As a consequence, the tensor ft,^^ 
has also two degrees of freedom. These degrees of freedom are related to clockwise and counter-clockwise rotation of 
matter. 



C. Scalar Perturbations 

Differentiation of equations (109b) covariantly with respect to we obtain 

^yfeiiij - ^N/cii ^ 2K(e(o) -l-p(o))u(i)||,|j. (137) 
Interchanging i and j in this equation, and subtracting the resulting equation from (137) we get 

'^tmj ~ ^lljIfcK " ^M^m +P(o))(W(i)p|j - U(i,||j|,), (138) 

where we have used that = ^lk\j\i- Using that V A = 0, we find from (111a) that the function C, must 

obey the equations 

C"K|fcb-C"b|feK = 0. (139) 

These equations are fulfilled identically in flrw universes. This can be seen as follows. We first rewrite these equations 
by interchanging the covariant derivatives in the form 

- C''\i\j\k) - (C'''|j|fc|i - C'^iMk) + (C'^'kij - C'''bii)|fc = 0- (140) 

Next, we use expression (117) and substitute the Riemann tensor (120) into the resulting expression. Using that 
CkIj = C\j\ii find that the left-hand sides of the equations (140) vanish. As a consequence, the equations (139) 
are identities. Therefore, the decomposition (111a) imposes no additional condition on the irreducible part of the 
perturbation k'j. 

We thus have shown that the system of first-order Einstein equations (109) is, for scalar perturbations, equivalent 
to the system 

Constraints: + ^ = -ke,!,, (141a) 

^k\^ - %|fe = + P(o))"(i)|K, (141b) 

Evolution: h\y + 3Hh\y + jHh\^ + 2 ^i^J^^n^. = -K,5'^{e^^ - p^,^), (141c) 
Conservation: + 3i7(£(i) + P(i)) + (e(o) + P(o))d(i) = 0, (141d) 

- 5(oV(i)|fc + 5ff(£(o) +Pm)ul^\\ = 0, (141e) 
+ iHn^r, + njojeid, = 0, (141f) 

where the local perturbations to the expansion, metric and the Ricci tensor are given by (73), (111a) and (114a) 
respectively. In the tensorial and vectorial case we found £(1, = and n^i^ = 0, implying that efl-^ = and nfl^ = 0, 
which made the tensorial and vectorial equations irrelevant for our purpose. Such a conclusion cannot be drawn from 
the equations (141). Perturbations with £(1, ^ and n^l•, ^ are usually referred to as scalar perturbations. 

Since the perturbation equations (141) contain only the components /ly^ , it follows that relativistic energy density 
and particle number density perturbations are characterized by two potentials, (j) 3.nd C- In the next section we will 
rewrite the system (141) into a form which is suitable to study the evolution of the quantities (3). 



i a r 

cdi. 



(£(0) +P(0))U(i) 



VIII. FIRST-ORDER EQUATIONS FOR SCALAR PERTURBATIONS 

The system of equations (141) can be further simplified by taking into account the decomposition (124). In the 
foregoing section we have shown that only the longitudinal part M(i)|| of the three- vector u^^ is coupled to density 
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perturbations. Using the properties (124)-(127) we can rewrite the scalar perturbation equations (141) in terms of 
quantities and which are suitable to describe exclusively the scalar perturbations. The result 

is that the first-order quantities occurring in the definitions (3) do explicitly occur in the set of equations. A second 
important result is that the metric components /ly^. occur, in the resulting equations, only in An evolution 

equation for this quantity follows from the (0, i) perturbed constraint equations and will be derived below. 

We now successively simplify all equations of the set (141) by replacing by its divergence (74), and 

eliminating /ip with the help of (73) and using that /I'^j. = as follows from (110) and (111). For equation (141a) 
we find 

2ff(0<,,-^<,,)-r^(i,ll ='^-£(1)- (142) 

Thus the (0, 0)-component of the constraint equations becomes an algebraic equation which relates the first-order 
quantities -iJ^), and e^iy In (142), is given by 

^^(i)!| = 9lo)ihllk\i\j ~ ^h\j\k) + 5 ^^(n)^'|/c' (143) 

as follows from (87) with (111)-(113) and (115). 

We will now derive an evolution equation for from equations (141b). Firstly, multiplying both sides of 

equations (141b) by gl^^ and taking the covariant divergence with respect to the index j we find 

.9(o)(^tfcHj - hlr\k\j) = 2K(e(o) +Pm)^w, (144) 



where we have also used (74). The left-hand side will turn up as a part of the time derivative of the curvature . In 

(0)ij ' 



fact, differentiating (143) with respect to ct and recalling that the connection coefficients Fj^^.j , (62), are independent 



of time, one gets 

^Ti-^^ii - -2H%,,^^ + 9%{hlm, - hUku) + 5 '^(0)%' (145) 
where we have used (28), (56), (90) and 

= (146) 

which is a consequence of = and the symmetry of hV . Next, combining equation (144) with (145), and, finally. 



ho)\k 

eliminating /iffj. with the help of (73), one arrives at 



%,,\\+2H%,,\\ -2«;(£(o, +p(0))i?(i, + i'i?,o,(e(i)-^9(i,) = 0. (147) 

In this way we managed to recast the three (0, z)-components of the constraint equations in the form of one ordinary 
differential equation for the local perturbation, to the spatial curvature. 

We now consider the dynamical equations (141c). Taking the trace of these equations and using (73) to eliminate 
the quantity /lyj,, we arrive at 

0,,, - + 6H{0,r> - ^?(i,) - %r,i\ = i«;(e<i, - P(i)). (148) 

Thus, for scalar perturbations, the three dynamical Einstein equations (141c) with i = j reduce to one ordinary 
differential equation for the difference 9^l^ — i^^). For i ^ j the dynamical Einstein equations are not coupled to scalar 
perturbations. Therefore, these equations need not be considered. 

Taking the covariant derivative of equation (141e) with respect to the metric g(o)ij and using (74), we get 

(£(0) +P(o))i9(i) -g(o)P(i)|fcK + 5i/(e(o)-(-P(o))^?(i) =0, (149) 



where we have used that the operations of taking the time derivative and the covariant derivative commute, since the 

^k 
(0)ij I 

form 



connection coefficients Fj^^.j , (62), are independent of time. With equation (92), we can rewrite equation (149) in the 
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where is the generalized Laplace operator which, for an arbitrary function f{t, x) and with respect to an arbitrary 
three-dimensional metric g^^{x), is defined by 

Thus, the three first-order momentum conservation laws (141e) reduce to one ordinary differential equation for the 
divergence 

Finally, the conservation laws (141d) and (141f) are already written in a suitable form. This concludes the derivation 
of the first-order equations for scalar perturbations. 

We thus have shown that the system of equations (141) is equivalent to the system 

Constraint: 2iJ(0(i) - i^d,) - \ = Ke^), (152a) 

Evolution: + 2H%,^\\ - 2Ke^„^{l + w)^^r, + | %„^{0^,J - ??(i,) = 0, (152b) 

9,,, - d,,, + 6H{e^,, - ,9,„) - = |K(e(„ - p,„), (152c) 

Conservation: £(i) +3H{e(i) + P(i)) +e(o)(l + w)6'(i) = 0, (152d) 

+ H{2 - 3/32)1?,,, + / = 0, (152e) 

+3i?n(i, +n(o,6l(i, = 0. (152f) 

The quantities /3(i) and wit) occurring in equations (152) are defined by 

m^^l^y Mt)^^. (153) 

The algebraic equation (152a) and the five ordinary differential equations (152b)-(152f), is a system of six equations 
for the five quantities 0^i^, ^R(i), £(1) and rt,,, respectively. This system is, however, not over-determined: in 
exactly the same way as we have shown in Section VI C 2 for the background equations, it can be easily shown by 
differentiation of the constraint equation (152a) with respect to time that the general solution of the system (152a)- 
(152b) and (152d)-(152f) is also a solution of the dynamical equation (152c). Therefore, it is not needed to consider 
equation (152c) anymore. 



IX. SUMMARY OF THE BASIC EQUATIONS 

In this section we will show that there exist unique gauge-invariant density perturbations (3). To that end we 
summarize the background and first-order equations for scalar perturbations. 



A. Zeroth-order Equations 

The Einstein equations and conservation laws for the background flrw universe are given by (57), (88), (90) and 
(92)-(93): 

Constraint: 3H'^ = ^ ^i?(o) + ^£(0) + A, (154a) 

Evolution: %„)^-2H%„), (154b) 

Conservation: £(„) = — 3iJ£(o)(l + w), (154c) 

^?(o, = 0, (154d) 

7=1(0) = -iHn^ny, (154e) 

where the initial value for (154b) is given by (91). As we have shown in Section VIC 2, the dynamical Einstein 
equation (89) is not needed. The set (154) consists of one algebraic and three differential equations with respect to 
time for the four unknown quantities e^o), '^■(o)! ^(o) = 3iJ and '^R^a)■ The pressure P(o)(0 is related to the energy 
density £(o)(i) and the particle number density n(of(t) via the equation of state (59). 
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We rewrite the Friedmann equation (154a) by dividing both sides by in the form 

1 = Ocurv + f^bar + ^cdm + ^lad + ^A, (155) 

where rJcurvi ^^ban f^cdmi ^lad E^^d Q\ are the contributions due to curvature; baryonic (ordinary) matter; CDM; 
radiation and dark energy respectively. The time dependence of these contributions is given by 

r!curv(i) = — ^^bar(^) + ^^cdm(^) + a•ad(^) = nA{t) = -^, (156) 

where we have used (65). This enables us to link the observations made with the wmap satellite to our treatise on 
density perturbations. The present day values of the quantities in (155) are, for a ACDM universe (also referred to 
as the concordance model) [14-18], given by 

^^bar(ip) = 0.0441, r!edm(tp) = 0.214, aad(ip)=0, ^^{t^) = Q.1A2. (157) 

With (155) we get 

f^curv(tp) = 0.000. (158) 

Thus, WMAP-observations indicate that the universe is flat. Moreover, it follows from WMAP-observations that the 
present value of the Hubble function (54) is 

-H(<p) = 71.9km/s/Mpc. (159) 

Using that 1 Mpc = 3.0857 x 10^^ m, we get for the curvature parameter k and the cosmological constant A, using the 
observed values (157)-(159), 

fc = 0, A = 1.34 X 10-^2 m-^ (160) 

respectively. The cosmological constant represents dark energy, also known as quintessence, a constant energy density 
filling space homogeneously. The existence of dark energy is postulated in order to explain recent observations that 
today the universe appears to be expanding at an accelerating rate. Since accelerated expansion takes place only at 
late times, we do not take into account A in our calculations of star formation in Sections XIII-XV. 

In astrophysics and cosmology, CDM is hypothetical matter of unknown composition that does not emit or reflect 
enough electromagnetic radiation to be observed directly, but whose presence can be inferred from gravitational effects 
on visible matter. One of the main candidates for CDM are the so-called wimps (Weakly Interacting Massive Particles) 
with a mass of approximately 10-10"^ GeV/c^. A recent estimate [38] yields a WIMP mass of approximately 70 GeV/c^. 
In comparison, the mass of a proton is 0.938 GeV/c^. Because of their large mass, WiMPs move relatively slow and 
are therefore cold. Since wimps do interact only through weak nuclear force with a range of approximately 10~^^m, 
they are dark and, as a consequence, invisible through electromagnetic observations. The only perceptible interaction 
with ordinary matter is through gravity. In the literature, CDM is therefore treated as 'dust,' i.e. a substance which 
interacts only through gravity with itself and ordinary matter. 



B. First-order Equations 

Since the evolution equation (152c) is not needed, as we have shown at the end of Section VIII, the first-order 
equations describing density perturbations are given by the set of one algebraic equation and four differential equa- 
tions (152a)-(152b) and (152d)-(152f): 

Constraint: 2H{9^i^ — ?9(i)) — \ ^R(i)\\ — (161a) 
Evolution: + 2H%,,\\ - 2Ke,„-,{l + w)i9(i) + | %„■,{0^^^ - ??(i,) = 0, (161b) 

Conservation: £(1, +3H{e(i) +e(o)(l + w)6'(i) = 0, (161c) 

^<,+i7(2-3/^^).,, + ^-^^=0, (161d) 

77(1) + 3i7n(i) + n(o)6'(i, = 0, (161e) 

for the five unknown functions e^j, n^i^, i?;!,, ^R(i-,\\ and 0^-^^ respectively. The first-order perturbation to the pressure 
is given by the perturbed equation of state (77) . 
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From their derivation it follows that the set of equations (161) is equivalent to the system of equations (141). The 
metric is in the set (161) contained in only one quantity, namely the local first-order perturbation '^R(i)\\ to the global 
spatial curvature. Using the sets of equations (154) and (161), we show in Section XII that our perturbation theory 
yields the Newtonian theory of gravity in the non-relativistic limit of an expanding universe with (3) as the key 
quantities. 



C. Unique Gauge-invariant Density Perturbations 



The background equations (154) and first-order equations (161) are now rewritten in such a form that we can draw 
an important conclusion. Firstly, we observe that equations (161) are the first-order counterparts of the background 
equations (154). Combined, these two sets of equations describe the background quantities and their corresponding 
first-order quantities: 

(£(0), £(!)), («(„,, n,i,), (^(0, = 0, t9(i,), ('^(0), (^(0) - 3i/, 0,,,). (162) 

Secondly, we remark that the quantities i3 and '^i? are not scalars of space-time, since their first-order perturbations 
"(9(1) and transform according to (C7) and (C9) respectively. In contrast, the first-order perturbations £(1), 

and 0(1) of the scalars (20) transform according to (15). Thus, only three independent scalars £, n and 6, (20), play a 
role in a density perturbation theory. Consequently, the only non-trivial gauge-invariant combinations which can be 
constructed from these scalars and their first-order perturbations are the combinations (3). We thus have shown that 
there exist unique gauge-invariant quantities £^) and nfly In Section XII we show that a perturbation theory based 
on these quantities yields the Newtonian theory of gravity in the non-relativistic limit of an expanding universe. 

By switching from the variables £(i), n^j, '^R(i)\\ and to the variables £^i) and nfl^ we will arrive at a set of 
equations for efl^ and nfl^ with a unique, i.e. gauge-invariant solution. This will be the subject of Section XI. First, 
we derive some auxiliary expressions related to the entropy, pressure and temperature. 



X. ENTROPY, PRESSURE, TEMPERATURE AND DIABATIC PERTURBATIONS 



In order to study the evolution of density perturbations, we need the laws of thermodynamics. In this section we 
rewrite the combined First and Second Laws of thermodynamics in terms of the gauge-invariant quantities (3) and 
the gauge-invariant entropy per particle sfly 



A. Gauge-invariant Entropy Perturbations 



Consider a gas of N particles with volume V. Let /j, be the thermodynamic — or chemical — potential, p the pressure 
and S its entropy. Then the internal energy E is given by the relation 

E = TS-pV + ^iN. (163) 

In terms of the energy per particle e = E/N, the entropy per particle s = S/N and the particle number density 
n = N/V this relation reads 

e^Ts-pn-^ + fi, (164) 

implying that 

de = Tds + sdT — n^^dp — pdn^^ + d/i. (165) 

From (163) we also find 

dE = TdS + SdT - Vdp - pdV + ^idN + Nd^. (166) 
The combined First and Second Laws of thermodynamics reads 

dE = TdS - pdV + AidTV. (167) 
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From (166) and (167) we find after division by N 

sdT ~ n"Mp + d/i = 0. (168) 

This relation enables us to eliminate d/i from (165). We so find 

Tds = de + pdn-^ . (169) 

This thermodynamic relation is independent of N and /i. In terms of the energy density defined as e = ne, we so find, 
finally, 

Tds = d( -) + pdf -) . (170) 



This is the relation we shall use in the following. 

The thermodynamic relation (170) is true for a system in thermodynamic equilibrium. For a non-equilibrium system 
that is 'not too far' from equilibrium, the equation (170) may be replaced by 

where d/dt is the time derivative in a local co-moving Lorentz system. Now, using e = £(o) -t- s = S(o) -I- S(i), 
P = P(o) +P(i) s-iid n = 71(0) + ^(1); we find from equation (171) 

where we neglected time derivatives of first-order quantities. With the help of equations (154c), (154e) and (153) we 
find that the right-hand side of equation (172) vanishes. Hence, S(o) — 0, implying that, in zeroth-order, the expansion 
takes place without generating entropy: s^o^ is constant in time. Hence, in view of (15), which is valid for any scalar, 
and S(o) = 0, the first-order perturbation S(i) is automatically a gauge-invariant quantity, i.e. S(i) = S(i), in contrast 
to £(1) and riji), which had to be redefined according to expressions (3). Apparently, the entropy per particle S(i) is 
such a combination of s^if and n^^ that it need not be redefined. This can be made explicit by noting that in the 
linear approximation we are considering, the combined First and Second Laws of thermodynamics (170) should hold 
for zeroth-order and first-order quantities separately. In particular, equation (170) implies 

T(o)S(i) = - ( £(1) ^'^(1) ) I (1'''3) 

"■(0) \ "-(0) / 

where we neglected products of differentials and first-order quantities, and where we replaced d£ and dn by £(i) and n^i^ 
respectively. We now note that the linear combination in the right-hand side of equation (173) has the property 

g(0) +-P(0) „ _pgi g(0) +P(0) si /1 74^ 



"(0) "(0) 

as may immediately be verified with the help of (3) and the equations (154c) and (154e). The right-hand side of 
expression (174) being gauge-invariant, the left-hand side must be gauge-invariant. This observation makes explicit 
the gauge-invariance of the first-order approximation to the entropy per particle, S(i). In order to stress the gauge- 
invariance of the correction S(i) to the (constant) entropy per particle, S(o), we will write sfl^, rather than S(i). From 
(173), (174) and sf'j = S(i) we then get 

'^(o)Sfi) = - (eu) ^"•fi) ) ' i^'^^) 

"(0) \ "(O) / 

where w is the quotient of zeroth-order pressure and zeroth-order energy density defined by (153). With (175) we 
have rewritten the thermodynamic law (167) in terms of the quantities efl^ and nfly 
We rewrite equation (175) in the form 

r,„,4 = -£i^iiil±i^^<), (176) 



where the gauge-invariant, entropy related quantity afl^ is given by 



£(0)(l + 7i;) 

The quantity afl^ occurs as the source term in the evolution equations (197) below. 



(177) 
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B. Gauge-invariant Pressure Perturbations 

We will now derive a gauge-invariant expression for the physical pressure perturbations. To that end, we first 
calculate the time derivative of the background pressure. From the equation of state (59) we have 

P(o) =^'n?^(o) +P6e(o), (178) 

where and p„ are the partial derivatives given by expressions (78) and (79). Multiplying both sides of this expression 
by 0(1) /^(o) and subtracting the result from (77) we get 

where we have used (3) to rewrite the right-hand side. Since p„ and depend on the background quantities 
and 71(0) only, the right-hand side is gauge-invariant. Hence, the quantity pf') defined by 

pfl^p,,,-P^0,,^, (180) 
is gauge-invariant. We thus obtain the gauge-invariant counterpart of (77) 

pfl^PeSfl+Pnnfl. (181) 

We will now rewrite this expression in a slightly different form. From (153) and (178) we obtain /3'^ — ~\~Pn (^(o) / "^(o) ) ■ 
Using equations (154c) and (154e) we find 

p'=p.+ (182) 

£(,-,)( 1 + W) 

With this expression and (177) and (182) we can rewrite the pressure perturbation (181) as 

pfl^f3hfl+pnafly (183) 

We thus have expressed the pressure perturbation pf^^ in terms of the energy density perturbation efl^ and the entropy 
related quantity af^^ rather than e^') and the particle number density perturbation nfly 

Expression (183) can be rewritten into an equivalent expression containing the entropy perturbation sfl^ explicitly. 
For pn we find 

^fdp\ ^fdp\ f(h\ £(o)(l + w) ^fdp 



dnj \dsj \dnj nf T^o) ' 



ys 



where we have used (175). Combining (176) and (183) we arrive at 

pfl=/3'efl+pssfl,. (185) 

Substituting (184) into (182), we get 

l3'=Ps-^^. (186) 

Expressions (182)-(183) refer to an equation of state p ~ p{n,e), whereas expressions (185)-(186) refer to the equiv- 
alent equation of state p = p(s, s). 

C. Gauge-invariant Temperature Perturbations 

Finally, we will derive an expression for the gauge-invariant temperature perturbation r(^) with the help of (A2a). 
For the time derivative of the energy density e(o)('^(o)i T^g^) we have 
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Replacing the infinitesimal quantities in (A2a) by perturbations, we find 



£(1) 



de 
df 



Multiplying both sides of (187) by 6'(i)/0(o) and subtracting the result from (188) we get 



t(l) 



de 
dn 



where we have used (3). Hence, the quantity 



de_\ 



^(1) — -'(I) A 



(188) 



(189) 



(190) 



is gauge-invariant. Thus, (189) can be written as 



dn J 



T 



9£\ 



(191) 



implying that T(^' can be interpreted as the gauge-invariant temperature perturbation. We thus have expressed the 
perturbation in the absolute temperature as a function of the perturbations in the energy density and particle number 
density for a given equation of state of the form e — s{n, T) and p — p{n, T). This expression will be used in Section XI 
to derive an expression for the fluctuations in the background temperature, St, a measurable quantity. 
Finally, we give the evolution equation for the background temperature T^ofit). From (187) it follows that 



-3iJ 



£(0,(1 -I- w) 



dn J J 



de_\ 



(192) 



where we have used equations (154c) and (154e). This equation will be used to follow the time development of the 
background temperature once £(o)(i) and f^(o)(^) are found from the zeroth-order Einstein equations. 



D. Diabatic Perturbations in a FLRW Universe 



In Section XA we have shown that the universe expands adiabatically. In this section we investigate under which 
conditions local density perturbations are adiabatic. 

By definition, an adiabatic, or isocaloric process is a thermodynamic process in which no heat is transferred to or 
from the working fiuid, i.e. it is a process for which 5Q — 0. For a reversible process we have 5Q = T(o)sf]'). Hence, 
reversible and adiabatic processes are characterized by TjajS^J) — 0. From expression (175) it follows that for adiabatic 
perturbations we have 

"(o)£a) - £(o)(l + w)nfl^ = 0. (193) 

Using the background conservation laws (154c) and (154e), we get the adiabatic condition for density perturbations 
in a FLRW universe 

ntojefj) - e^o^n% = 0. (194) 

In a non-static universe we have e^o) 7^ and n(o) ^ 0. In this case, equation (194) is fulfilled if and only if the energy 
density is a function of the particle number density only, i.e. if e = e{n). In particular, density perturbations in a 
perfect pressureless fluid with e = nmc^ are adiabatic. This is the case in the non-relativistic limit w/c of an 
expanding FLRW universe, see Section XII. In all other cases, e = e{n, T) and p = p{n, T), local density perturbations 
evolve diabatically. 
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XI. MANIFESTLY GAUGE-INVARIANT FIRST-ORDER EQUATIONS 

The five perturbation equations (161) form a set of five equations for tlie five unknown quantities S(i^, 

and ^^ly This system of equations can be further reduced in the following way. As has been explained in 
Section IV B, our perturbation theory yields automatically 9fl^ — 0, (25). As a consequence, we do not need the 
gauge dependent quantity 6(iy Eliminating the quantity ^(i) from equations (161) with the help of equation (161a), 
we arrive at the set of four first-order differential equations 

1 f 



£(1) + 3i7(£(i) + P(i)) + £(0,(1 + w) 

71(1) + 3i?n(i) + n(o) 



2H 
1 

2H 



1 3 
2 



( K£(i) 



1 3n 

2 ^(1) 



1 



^-+^(2-^^^)^--^£,,(l + .) 

- 2k£(o,(1 + U;)7?(i, + (k£(i, + i 



= 0, 
= 0, 

= 0, 



(195a) 
(195b) 
(195c) 
(195d) 



for the four quantities £(i), rijij, iJjij and From their derivation it follows that the system of equations (195) is, 

for scalar perturbations, equivalent to the full set (109) of first-order Einstein equations. 

The system of equations (195) is now cast in a suitable form to arrive at a system of manifestly gauge-invariant 
equations for the physical quantities £^) and nfl^ since we then can immediately calculate these quantities. Indeed, 
eliminating the quantity 6^i^ from equations (3) with the help of equation (161a), and using the background equations 
(154) to eliminate the time derivatives £(o), n^g^ and 9^o) = 3H, we get 



^i?(o) + 3k£(o)(1 + w) 

3n(o)(K£(i) + 2^7^(1) + F^(i) 
^R^o) + 3k£(o)(1 + w) 



(196a) 
(196b) 



The quantities sfl^ and nfl^ are now completely determined by the system of background equations (154) and the 
first-order equations (195). 



A. Evolution Equations for Density Perturbations 



Instead of calculating sfl^ and nfl^ in the way described above, we proceed by first making explicit the gauge- 
invariance of the theory. To that end, we rewrite the system of four differential equations (195) for the gauge 
dependent variables £(i), riji), li^i^ and into a system of equations for the gauge-invariant variables £^') and nfly 

It is, however, of convenience to use the entropy related perturbation <jfly defined by (177), rather than the particle 
number density perturbation nfly The result is 



- ai£(i) -t- a2£(i) — 03(7(1), 



(197a) 

The derivation of these equations is given in detail in Appendices B 1 and B 2. The coefficients ai, a2 and 03 occurring 
in equations (197) are given by 



02 = AC£(o,(l + u.) - 4i7| + 2H^{2 - 3/32) + ^R^^ 



^R(o) + 3k£(o)(1 + w)J' 

2 



^ 5ff2(i + 3^2)_2i7^ 



■^i?(0) + 3k£(o)(1 + w) 



\ 



0-3 



-18i?2 



'%o) + 3k£(o)(1 + w) 



e(0)Pe«(l + w) + - (3'^Pn+PePn + n(o)P„„ 

oil p 




1 3n 

2 ^(0) 



1 3n 



(198a) 
(198b) 
(198c) 
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where the functions f3{t) and w{t) are given by (153). In the derivation of the above results, we used equations (154). 
The abbreviations p„ and are given by (78). Furthermore, we used the abbreviations 



Pn 



_cPp_ 
de dn 



(199) 



The equations (197) contain only gauge-invariant quantities and the coefficients are scalar functions. Thus, these 
equations are manifestly gauge-invariant. In contrast, the equations (195), being linear Einstein equations and con- 
servation laws, are themselves gauge-invariant, but their solutions are gauge dependent, see Appendix C for a detailed 
explanation. 

The equations (197) are equivalent to one equation of the third-order, whereas the four first-order equations (195) 
are equivalent to one equation of the fourth-order. This observation reflects the fact that the solutions of the four 
first-order equations (195) are gauge dependent, while the solutions efi, and afl^ of equations (197) are gauge-invariant. 
One 'degree of freedom,' say, the gauge function tp{x) has disappeared from the scene completely. 

The equations (197) constitute the main result of this article. In view of (177), they essentially are two differential 
equations for the perturbations eflj and nfl^ to the energy density £(o) (t) and the particle number density n(o) (t) 
respectively, for flrw universes with k = —1, 0, -1-1. They describe the evolution of the energy density perturbation ef'j 
and the particle number density perturbation nfl^ for FLRW universes filled with a fiuid which is described by an 
equation of state of the form p = p{n, e), the precise form of which is left unspecified. 



B. Evolution Equations for Contrast Functions 



In the study of the evolution of density perturbations it is of convenience to use a quantity which measures the 
perturbation to the density relative to the background densities. To that end we define the gauge-invariant contrast 
functions 5^ and (5„ by 



Using these quantities, equations (197) can be rewritten as (Appendix B3) 

5. 



1 + w 



1 d 



iHn^o)Pn 



1 + w 



c dt \ 1 + w J e(o)(l + w) 
where the coefficients bi, 62 and 63 are given by 

H p \6H ■^i?(o) -I- 3Ke(n)(l + w) 

62 = -iKe(o)(l + w){l + 3w) + H'^ {l-3w + 6l3^{2 + 3w)) 
Ke,o)(l + w) ^ ^ ""2 



( w + TT , ^ 

P\ 3^(0, -H3k£(o)(1-Hw) 



(0) 2 



^i?(o, + 3k£(o)(1 + w) J \a 



63 = 



%(0) + 3k£(o)(1 ■ 



£(0)P£n(l +'^> + P Pn+PePn + ?^(o)^lr, 



+ Pn 



(200) 

(201a) 
(201b) 

(202a) 




In Sections XIII-XV we use the equations (201) to study the evolution of small energy density perturbations and 
particle number perturbations in flrw universes. 

The combined First and Second Laws of thermodynamics (175) reads, in terms of the contrast functions (200), 



T s^' 

-t(0)*(l) 



e(o)(l + w) f _ 5^ 
" 1 + w 



(203) 



where sf\ is the entropy per particle. 
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Finally, we define the relative temperature perturbation 5t and the relative pressure perturbation 5p by 



Srit,x)^^^^, S^it^x)^-^^^^. (204) 
Using the expressions (191), (200) and (204) we arrive at the relative temperature perturbation 

(205) 



de_\ 



The relative pressure perturbation follows directly from (181). We get 

£(0) f dp\ J , 71(0) f dp 



= — + — M- (206) 

We thus have found the relative temperature and the relative pressure perturbations as functions of the relative 
perturbations in the energy density and particle number density for an equation of state of the form e = e{n,T) and 
p = p{n, T) (see Appendix A). The quantity (5t(^, x) is a measurable quantity in the cosmic background radiation. 



XII. NON-RELATIVISTIC LIMIT IN AN EXPANDING FLRW UNIVERSE 

In Section IX C we have shown that there exist only two gauge-invariant quantities e^^, and n^l^ , which could be 
the real energy density and particle number density perturbations. In this section we show that in the non-relativistic 
limit v/c — ^ the quantities e^^) and nfl^ become equal to their Newtonian counterparts. This implies that e^^ and 

nfl^ are indeed the local perturbations to the energy density and particle number density perturbations. 

It is well known that if the gravitational field is weak and velocities are small with respect to the velocity of light 
(v/c — !■ 0), the system of Einstein equations and conservation laws may reduce to the single field equation of the 
Newtonian theory of gravity, namely the Poisson equation (229). In the first-order perturbation theory developed in 
this article, the gravitational field is already weak, so that, at first sight, all we have to do to arrive at the Newtonian 
theory of gravity is to take the non-relativistic limit w/c — > in all equations. Since in the Newtonian theory the 
gravitational field is described by only one, time-independent, potential fix), one cannot obtain the Newtonian theory 
by simply taking the non-relativistic limit v/c ^ 0, since in a relativistic theory the gravitational field is described, 
in general, by six potentials, namely the six components hij{t,x) of the metric. In this article we have used the 
decomposition (110). Moreover, we have shown in Section VII that only given by (111a) is coupled to density 
perturbations. By using this decomposition we have reduced the number of potentials to two, namely (j)(t,x) and 
(^{t,x). We have rewritten the system of equations (141) for scalar perturbations into an equivalent system (161). As 
a result, the perturbation to the metric, hi-, enters the system (161) via the trace 



^(1)1! ^2 



\\3 

2 



~ Cl'l'i^ife + C"|J'|; + I %o>i^^ + C"|fc)l , (207) 



of the spatial part of the perturbation to the Ricci tensor (114a), and via the perturbed expansion scalar (73) 

0a)=^(i,-^(3(/' + C'V-)' (208) 

where we have used (Ilia). This shows explicitly that density perturbations are, in flrw universes, described by two 
potentials 0(t, a;) and ({t,x). In this section we show that in the non-relativistic limit of a flat flrw universe, the 
potential ( drops from the perturbation theory. For a fiat (i.e. fc = 0) flrw universe we have '^R(o) = 0, as follows 
from (65). The perturbation to the spatial part of the Ricci scalar, (207), reduces in this case to 

4 ,,, 4 V^(h , , 

'^u,||-^</'"|fc = -^^, (209) 

where is the usual Laplace operator see (51) and (151). For a fiat flrw universe, the zeroth-order equations (154) 
reduce to 

Constraint: 3iJ^ = Ke^a) + A, (210a) 
Conservation: = — 3i/e(o)(l -l- w), (210b) 

ri(o) = -3Hn^o)- (210c) 
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Upon substituting (209) into (161) and putting ^R^o) = 0, we arrive at the set of perturbation equations for a flat 
FLRW universe: 

1 V^(t) AttG 



Constraint: H{e,,, ^ d,,,) + ^ ^[e^l, + f^6,,A, (211a) 



Evolution: ^ + ^£(0,(1 + w)dn. = 0, (211b) 

Conservation: £(1, + 3if(e(i) + P(i)) + £(0) (1 + w)^(i) = 0, (211c) 

+ H{2 - 3/?2)^,„ + I = 0, (211d) 

+ 3i?n(i) + 7i(o)6'(i) = 0, (211e) 

where we have used (38), and our new definition (3a) to ehminate e^l) from (161a). The scale factor of the universe 
a{t) follows from the Einstein equations (210) via H = a/a. In the equations (211), the potential C, occurs only in the 
quantity 9^^, see (208). 

We now consider the sets of equations (210) and (211) in the non-relativistic limit v/c^ 0. Since the spatial part 
UjjjII of the fluid four-velocity is gauge dependent, (C6b), with a physical component and a non-physical gauge part, 
we deflne the non-relativistic limit u/c — > by 

W „ 

I ^ (1)11 physical , . 

"(i)ilphysical — ^ ^' l^-L^J 

i.e. the physical part of the spatial fluid velocity is negligible with respect to the speed of light. In this limit the 
kinetic energy per particle ^m{v'^) = ^kBT — is small compared to the rest energy mc^ per particle, implying that 
the pressure p = nk^T — ^ (n 7^ 0) vanishes also and that the energy density of the universe is given by £ = nmc^. 
In other words, the non-relativistic limit (212) implies 

p = 0, e ~ nm(? , w=-=0. (213) 

£ 

In view of (213), the background equations (210) take on the simple form in the non-relativistic limit (212) 

Constraint: 3H^ — ne^a) + A, (214a) 
Conservation: £(„) — --3i/£(o), (214b) 
h^o) = -SHn^oy (214c) 

Thus, even in the non-relativistic limit, a non-empty universe cannot be static, therefore we have H ^ 0. Note that 
in the limit (212), equations (214b) and (214c) are identical, since then £(o) — n^g^mc^, in view of (213). 
The consequences of the limit (212) are as follows. In the limit p — 0, equation (141e) reads 

^^(£(o)Wa)ll) =0- (215) 

Using the background equation (214b), equation (215) takes on the simple form 

= -2i/,i;,,ii. (216) 

The general solution of this equation is, using H = a/a, 

<rM.u,.--^f'd,^i^). (217) 

Since the physical part of u^^^y vanishes in the limit (212), the solution (217) is a gauge mode, as follows from (C6b). 
Thus, in the limit (212) we are left with the non-physical quantity (217). If a quantity is proportional to the gauge 
function ip{x) or its partial derivatives, it may be put equal to zero without loss of physical information. If we require 
that in the limit (212) U(i)||gaugo ~ ^ holds true before and after a gauge transformation we find from (217) that 
dk'4>{x) = or, equivalently. 



V'(a;) = V'- (218) 
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In view of (218), the gauge dependent functions and n^iy transform under a gauge transformation in the limit 
(212) according to (2), with constant ip: 

£(1) ^ £(1) = £(1) + ^'£(0), (219a) 
71(1, ^- 71(1, = n,!) + -071(0). (219b) 

Since in an expanding universe the time derivatives £(o)(^) 7^ 0, (214b), and 7i(o,(t) 0, (214c), the functions £(i,(i:,x) 
and 71(1, (t, a;) do still depend upon the gauge, so that these quantities have, also in the limit (212), no physical 
significance. Finally, it follows from (C2) and (218) that the gauge transformation (4) reduces in the limit (212) to 
the gauge transformation 

x"~,x"~tp, x' -^x'-x\x). (220) 

In other words, time and space transformations are decoupled: time coordinates may be shifted, whereas spatial 
coordinates may be chosen arbitrarily. Thus, in the limit (212) the general relativistic infinitesimal coordinate trans- 
formation (4), with given by (C2), reduces to the most general infinitesimal coordinate transformation (220) which 
is possible in the Newtonian theory of gravity. By now, it should be clear that the gauge problem of cosmology cannot 
be solved by 'fixing the gauge,' since the constant tp and the three functions x*(^) cannot be determined. 

We now consider the perturbation equations (211) in the limit (212). Upon substituting i?(i, = (ttj^^j||)|fc = 0, 
w = P(o)/£(o) = and P(i) = into these equations we find that equation (211d) is identically satisfied, whereas the 
remaining equations (211) reduce to 

Evolution: V^^ = 0, (221b) 
Conservation: e(i, + 3i?e(i) + £(o)^(i) = 0, (221c) 
71(1, + 3i/77(i, + n(o,6'(i) = 0. (221d) 

The constraint equation (221a) can be found by subtracting ^6^i-,/H times the time-derivative of the background 
constraint equation (214a) from the constraint equation (211a) and using that O^q^ — 3H, (57). This shows explicitly 
that (3a) is the only possible choice to construct the gauge-invariant quantity efly Equations (221c) and (221d) 
are identical in the non-relativistic limit (212) since £(i, = n^i^mc^, in view of (213). The quantities £(1, and 77(1,, 
which are gauge dependent in the General Theory of Relativity, are also gauge dependent in the non-relativistic limit 
(212), so that equations (221c) and (221d) need not be considered: the equations (221c) and (221d) have no physical 
significance and can be discarded. Consequently, the perturbed expansion scalar f?(i, does not occur anymore in the 
perturbation theory, and we are left with one potential only. 
Equations (221a) and (221b) can be combined to give 



Constraint: V^^ = a-^Ef' (221a) 



\/'4,{x)^^a'{t)efl{t,x). (222) 



Or, equivalently. 



VV(a;) = ^a'(ip)£f;,(tp,=r), (223) 



where tp indicates the present time. This Einstein equation can be rewritten in a form which is equivalent to the 
Poisson equation, by introducing the potential ip 



^(x) - (224) 



Inserting (224) into (223) we obtain the result 



VV(3^)=4^G ^'''^V^^ - (225) 

Finally, we have to check the expressions (3) or, equivalently, (196) in the limit (212). It can easily be verified that 
in the non-relativistic limit (212) of a flat flrw universe, expression (196a) reduces to 



1 



4) = -^Xi,l|- (226) 
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In view of (38) and (209), this equation is equivalent to (221a). Using (226), we find that expression (196b) reduces 
in the non-relativistic Hmit (212) to 



'(1) — "(1) 



^(%,-4)- (227) 



Combining £(„) = n^o^mc^ and £(i) = n^jmc^, (213), with (227) yields 

= (228) 

which is the gauge-invariant counterpart of e^iy — n^-i^simc^ . Again, we have to conclude that (3b) is the only possible 
choice to construct the gauge-invariant quantity nf^y 

Combining equations (225) and (228), we arrive at the Poisson equation of the Newtonian theory of gravity: 

VV(=c) = 47rG/9(i)(a;), (229) 

where 

p,„(a;)^?^^^^^, p,,,{x)^mni,{t,,,x), (230) 

is the mass density of a perturbation. 

In this section we have shown that in the non-relativistic limit (212) of a flat FLRW universe, the first-order 
perturbation equations (161) together with the new definitions (3) reduce to the well-known Newtonian results (220), 
(228) and (229). Consequently, the gauge-invariant quantities ef^ and n^'j given by (3) are indeed the energy density 
and particle number density perturbations. 



XIII. PERTURBATION EQUATIONS FOR A FLAT FLRW UNIVERSE 

In this section we derive the perturbation equations for a fiat flrw universe with a vanishing cosmological constant 
in its radiation-dominated, plasma-dominated and matter-dominated stages. 

The six numerical values that we need in the following are: i. the redshift at time toq when the matter density has 
become equal to the radiation density, ii. the redshift at time idee when matter and radiation decouple, in. the present 
value of the Hubble function, iv. the present value of the background radiation temperature, v. the age of the universe 
and, vi. the fiuctuations in the background radiation at decoupling. They follow from WMAP [14-18] observations: 

z(toq) = 3176, (231a) 

z(idcc) = 1091, (231b) 

H(tp) = 71.9km/sec/Mpc = 2.33 x 10-i*^s-\ (231c) 

r(o)^(tp) = 2.725 K, (231d) 

^ 13.7 Gyr = 4.32 x lO" sec, (231e) 

^T,(^dcc)<10-^ (231f) 

where we have used that 1 Mpc = 3.0857 x 10^^ m (1 pc = 3.2616 ly). 
The cosmological redshift z{t) is given by 



z{t) = ^~l, a(tp) = l, (232) 
where we have normalized the scale factor a{t) to unity at i = ip, which is only allowed in a flat flrw universe. 



A. Radiation-dominated Phase 



We consider the universe after the era of infiation. Let tj-ad be the moment at which exponential expansion of the 
universe has come to an end, and the radiation-dominated era sets in. The radiation-dominated universe starts with 
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a temperature T(o)7(tiad) ~ lO'^^ K and lasts until matter-energy equality, which occurs at the redshift z{tcq) = 3176. 
If we assume [see Weinberg [35], (15.5.7) for an explanation] that during the time that matter and radiation were in 
thermal contact the temperature of the radiation, r^o)^, dropped according to the formula T(o)^(i) — A/a(t), where A 
is a constant and a{t) the scale factor, we have 

^(0)7 (*) _ «(ip) (233) 



T(o)j{tp) a{t) ■ 

This relation follows also from the Einstein equations in Section XIII A 1 . We suppose that matter and radiation are 
still in thermal equilibrium at the end of the radiation-dominated era, i.e. at a redshift of z(teq) = 3176, yielding 
a{tp)/a{tcq) = 3177. In this way we find from (231d) and (233) 

T<„)^(teq) = 8657K. (234) 

Hence, as a rough estimate we may conclude that the universe is dominated by radiation in the temperature interval 

10^^ K > r(o)^ > 10"* K. (235) 

If we neglect, in the radiation-dominated era, the contribution of the electrons and neutrinos, we have for the photon 
energy density e 



06^(0)^, (236) 



where ae is the black body constant 



flB = = 7.5658 X 10"^^ J m-^ K^^ (237) 

15n.'^c'' 

with fee and h Boltzmann's and Planck's constant respectively. Furthermore, we could neglect the cosmologi- 
cal constant A, (160), with respect to the energy content of the universe at the matter-radiation equality time 
KaBT|o)^(icq) = 1-0 X 10~'*^m~^, see equation (154a). Finally, we consider a flat universe {k = 0), implying, with 
(65), that %o){t) = 0. 

1. Zeroth- order Equations 

In the radiation-dominated era, the contributions to the total energy density due to baryons, nm^c^, can be 
neglected, so that the equations of state read 

e(n, T) = a^T^, p{n, T) = iasT^. (238) 

Moreover, we can put A = for reasons mentioned at the end of Section IX A. The zeroth-order equations (154) then 
reduce to 

ff2 = iK£(o), (239a) 

£(0, - -4H£(o), (239b) 

71(0) = —^Hn^oy (239c) 



These equations can easily be solved, and we get the well-known results 



H{t) = \ (ct)-^ = iJ(t,ad) (t^) , (240a) 



£(«,(<) = |:(rf) ' (240b) 

«<o,(t)=n(o,(irad) f^^) . (240c) 

The initial values -ff(trad) and £(o)(irad) are related by the constraint equation (239a) taken at <rad- the time at which 
the radiation-dominated era sets in. 
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Using the definition of the Hubble function H = a/a we find from (240a) that 

a{t) = a{t,,d) (t^) ' . (241) 

V ^rad / 

Combining expressions (236), (240b) and (241), we arrive at (233), which can be rewritten in the form 

T(o,^(i) = r<o„(tp)[z(0 + 1], (242) 

where we have used (232). 

With (240)-(241) the coefficients bi, &2 and 63, defined by (202), occurring in the equations (201) can be calculated. 
This will be the subject of the next section. 



2. First-order Equations 

In view of (78), we find from (238) 

Pn=0, p, = l, (243) 

so that, according to (153), we have w = ^ and /3 = 1/a/3, see (182). The zeroth-order solutions (240) can now be 
substituted into the coefficients (202) of the equations (201). Since V = V for a ffat universe, the coefficients bi, 62 
and 63 reduce to 



bi = -H, 62 = --^ + |k£(„„ &3 = 0, (244) 



1 V2 

3 a 

where we have used (239a). For the first-order equations (201) this yields the simple forms 

Se~HS,-(^~- §^£(0)) 5, = 0, (245a) 

~{Sr.-l6,)= 0, (245b) 

where H = d/a, (55). Since the right-hand side of (245a) vanishes, the evolution of density perturbations is inde- 
pendent of the actual value of (5„ — '^S^, see (201) and (203). In other words, the evolution of density perturbations 
is not affected by perturbations in the entropy. Equation (245b) implies that (5„ — |(5e is constant in time. Conse- 
quently, during the radiation-dominated era, perturbations in the particle number density are gravitationally coupled 
to radiation perturbations: 

6n{t, X) - |(5e(t, X) = (5„(irad, x) - |(5e (trad , £c) , (246) 

where the right-hand side is constant with respect to time. From (203) we find for the perturbation of the entropy 
per particle: 

Sfl{t,x) = 4(irad,a;) = -^^%^k(irad,a)) " |<Se(trad)l, (247) 

'^(0) V''i'ad j L 

where we have used that TJo)^ cx a~^, n^o) oc a~^ and (246). Thus, throughout the radiation-dominated era the 
entropy per particle is constant with respect to time. Since n^,-,^ ^ and nfl^ ^ 0, we have sfl^ ^ 0, see (194). 
Equation (245a) may be solved by Fourier expansion of the function 6^. Writing 

S,{t,x)=5,it,q)e"i--, (248) 

with q — \q\ — 2tt/X, where A is the wavelength of the perturbation and i^ — —1, we find 

V^S,{t,x)^-qH,{t,q), (249) 

so that the evolution equation (245a) for the amplitude S^{t,q) reads 
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where we have used (239a)-(240). This equation will be rewritten in such a way that the coefficients become dimen- 
sionless. To that end a dimensionless time variable is introduced, defined by 



t 

trad 



t > t 



rad : 



with irad the time immediately after the inflationary era. This deflnition implies 



1 



C^iad 



" d" d" 
dr" ^ ^ dr" 



71=1,2,. 



(251) 



(252) 



where we have used (240a). Using (240a), (251) and (252), equation (250) for the density contrast 6^{T,q) can be 
written as 



s': 



1 

2^^ 



Z 1 

Mr 1 

4t 2t2 



where a prime denotes differentiation with respect to t. The constant /ir is given by 

q 1 1 

fJ-T 



(253) 



(254) 



a(trad) H{t,^d) \/3 ' 

The general solution of equation (253) is a linear combination of the functions J_|_ i (/ii.^/T)T'^/'*, where J^i{x) = 
y/2/{Trx) since and J_i{x) = y/2/{'Kx) cos x are Bessel functions of the flrst kind: 



Se{T, q) = Ai{q) sin (^t.%/r) + ^2(9) cos (/ir\/T) %/t, 
where the functions Ai{q) and j42(q) are given by 



Ai{q) = 4(irad,g)sin^i. 



^2(9) = 5e(trad,g)cOS/ii. 



COS fJL-c 
Air 

sin /ij. 

Mr 



'^£(^rad,g) - 
'^£(irad, g) - 



'^£(<rad,g) 
^f(irad) 

'^^(^rad.g) 
^(^rad) 



where we have used that 



^e(trad, ?) = ^^(r = 1, q), <5e(trad, = 2i/ (t,.ad )<5; (t - 1, g). 



(255) 

(256a) 
(256b) 

(257) 



as follows from (252). The relative perturbations in the particle number density, (5„, evolve by virtue of (246) also 
according to (255). 

We consider the contribution of the terms of the Fourier expansion of the energy density perturbation [see (248)] 
in two limiting cases, namely the case of small A (large q) and the case of large A (small q). 

For large-scale perturbations, A — > oo, the magnitude of the wave vector |q| = 27r/A vanishes. Writing 
5e{t) = 5e{t, q = 0) and 4(0 = 4(i, Q = 0), we flnd from (254)-(256) that, for t > U^d, 



H{tia,d) 



t 

^rad 



2<5e(irad) 



ad) 



H{trad) 



t 

^rad 



(258) 



It is seen that the energy density contrast has two contributions to the growth rate, one proportional to t and one 
proportional to t^/^. These have been found by a large number of authors. See Lifshitz and Khalatnikov [6], (8.11), 
Adams and Canuto [7], (4.5b), Olson [8], page 329, Peebles [9], (86.20), Kolb and Turner [10], (9.121) and Press and 
Vishniac [11], (33). The precise factors of proportionality, however, have not been published earlier. From the first 
of them we may conclude, in particular, that large-scale perturbations only grow if the initial growth rate is large 
enough, i.e. 



4(irad) > 5,{trad)H{t,.,A) ^ 5'^{t = 1, q) > \5,{t = 1, q). 



(259) 



otherwise the perturbations are decaying. For CDM and ordinary matter the same growth rate 5„ oc t is found in the 
literature for super-horizon perturbations, see, for example, the textbook of Padmanabhan ]39], Section 4.4. Thus, 
for large-scale perturbations, our treatise corroborates the outcomes found in the literature on the subject. 
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We now come to the second case. In the small-scale limit A 
(256), that 



(or, equivalently, \q\ oo) we find, using (254)- 



t 

^rad 



Mr - Mr 



t 

^rad 



(260) 



We see that in the limit of small A, the contribution to the growth rate is smaller than the leading term in the 
expression (258). Physically, this can be understood: on small scales, the pressure gradients |Vp| ~ p/X are much 
higher than on large scales. 

In Section XVI we review the standard results. Comparing our result (260) with the standard result (352) found 
in the literature, we observe that we obtain oscillating solutions with an increasing amplitude, whereas the standard 
equation (350) yields oscillating solutions (352) with a decreasing amplitude. 

The above calculated behavior of density perturbations in the radiation-dominated universe is important for star 
formation in the era after decoupling of matter and radiation: the oscillating growth shows up in the cosmic background 
radiation as random fluctuations on different scales and amplitudes (i.e. intensities in the background radiation). As 
we will show in Section XV, density perturbations yield massive stars. 

B. Plasma Era 

The so-called plasma era sets in at time t^q, when the energy density of ordinary matter equals the energy density of 
radiation, i.e, when n^o){teq)m'C^ = aBT^oni^m) ^ ^^'-^ ends at time tdoc, the time of decoupling of matter and radiation. 
In the plasma era the matter-radiation mixture can be characterized by the equations of state, see Kodama and 
Sasaki [25] Chapter V, 



s{n, T) = nmc^ + a^T^, p{n, T) = ^aBT^, 



(261) 



where we have not taken into account the contributions to the pressure from ordinary matter, m = mn, or CDM, 
m = mcdm, since these contributions are vanishingiy small in comparison to the radiation energy density. Using (A5), 
we find 



Pn 



Pe 



(262) 



Substituting the expressions (261) and (262) into the entropy equation (201b) and integrating the resulting equation 
yields 



l + w{t) 



^6 (^cqi x) 
1 + W{tcq) 





l*ct 


exp 









dr 



(263) 



It is well-known that ordinary matter perturbations are coupled to perturbations in the radiation density. This 
coupling is attributed to the high radiation pressure in the radiation-dominated and plasma eras. In addition to this 
coupling, expression (263) shows that perturbations in the particle number density, S„, are also gravitationally coupled 
to perturbations in the total energy density, S^- This has considerable consequences for the growth of perturbations in 
CDM. Since the energy density of CDM can be written as n(o,mcdmC^i (263) implies that also density perturbations in 
CDM are gravitationally coupled to perturbations in the radiation energy density. Consequently, CDM perturbations 
can, just as perturbations in ordinary matter, start to grow only after decoupling. In other words, CDM and ordinary 
matter behave gravitationally in exactly the same way. This may rule out CDM as a means to facilitate the formation 
of structure in the universe. The same conclusion, on different grounds, has also been reached by Nieuwenhuizen 
et al. [5]. 

C. Era after Decoupling of Matter and Radiation 

In the era after decoupling of matter and radiation, we have to distinguish between the matter temperature and 
the radiation temperature. The radiation temperature T(o)-y evolves as (233), whereas the matter temperature T(o) 
evolves according to (277). Once protons and electrons recombine to yield hydrogen, the radiation pressure becomes 
negligible, and the equations of state reduce to those of a non-relativistic monatomic perfect gas [Weinberg [35[, 
equations (15.8.20) and (15.8.21)[ 



(264) 
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where fee = 1.3806504 x 10"^'^ JK~^ is Boltzmann's constant, mn the mass of a proton, and T the temperature of the 
matter. Since the energy density in (264) is not of the form e = s{n), see Section XD, density perturbations cannot 
be adiabatic. 



1. Zeroth-order Equations 



The maximum gas temperature occurs around time tdcc of the decoupHng of matter and radiation and is equal to 
the radiation temperature T'(o)-y(idcc)- Using (231)-(233), we get for the temperature at the time of decoupUng 

r(„)(tdcc) - T(„,^(tdoc) - 2976 K. (265) 
Since the universe cools down during its expansion, it follows from (264) that 



P ^ MkW < Mkl^^iS^ = 2.73 X 10-10, t>idec. (266) 
e tohc^ TTifjc^ 

Hence, the pressure is negligible with respect to the energy density. This implies that, to a good approximation, 
£(0) ±J5(o) ~ £(0) Ei-nd £(0) ~ n(o)TOHC^. Hence, in an unperturbed flat flrw universe the pressure can, in the background 
equations, be neglected with respect to the energy density. The above facts yield that the Einstein equations and 
conservation laws (154) for a flat flrw universe reduce to 

i/2^iKe(o), (267a) 
£(0) = -3i/£(o), (267b) 
71(0) = -3i/n(o), (267c) 

where we have put the cosmological constant A equal to zero. The solutions of (267) are 

Hit) = I (ct)-^ = i/(imat) (j^) , (268a) 



^ ^mat / 

^<o)(i) = ^ (ci)"' = e(o,(imat) (t^) , (268b) 

-3 



a(imat) 



"-(0)(i) = "-(0)(^mat) ( V ) ■ (268c) 



The scale factor a{t) after decoupling, however, has a time dependence which differs from that of the radiation- 
dominated era (241). From (268a) we find 



2 

t \ 3 



a{t) = a(W) [- , (269) 

V Emat / 

where tmat is some initial time after decoupling of matter and radiation: 

idee < ^mat < ^p. (270) 

The initial values -ff(tmat) and £(o)(tmat) are related by the constraint equation (267a) taken at t — tmat. With 
(268)-(269) the coefficients (202) of the perturbation equations (201) are known functions of time. 



2. First-order Equations 



We first remark that, in the study of the evolution of density perturbations, we may not neglect the pressure with 
respect to the energy density. The case of a pressureless perfect fluid is already thoroughly discussed in Section XII on 
the non-relativistic limit. We neglect k^T^g^ / {m^c^) with respect to terms of the order unity in the final expressions. 

Using equations (A5) we find from the equations of state (264) 

Pe = l, p„ = -|mHC^ (271) 



From (182) it follows that 

m 



c2 + |A:Br(o,(t) 



tohc 



Since kBT^(,y(t) <C m]iC^, (266), we find, to a good approximation 



c y 3 tohc^ 

where Vs is the speed of sound. Differentiating (273) with respect to time yields 

For the time development of the matter temperature it is found from (192) and (264) that 
Combining (274) and (275) results in 



For the evolution of the matter temperature it is found from (275) that 



-^(0)(0 — ^(0)(^mat) \ ^fj^ \] — T'{0)(^mat 



where we have used that H = a/a and (232). 

We now consider equation (201b). Using (271) it is found for this equation 

^^(5n-'5s) = -2ff(5„-,5,). 
c at 

The general solution of equation (278) is, using also (269), 

A(t, x) = A(i„,at, x) f ^5^) = A(t^at, a;) (-^ 

\0.\tmat) J \Emat 

where the quantity A(t, a;) is defined by 

A(t, a;) EE (5„(i,a;) - 5^(t,x). 
The combined First and Second Laws of thermodynamics (203) reads for an equation of state (264) 

T{o)sfl) = -[rtiuc^ + |fcBT((,))A - fcBT'(o)<5„. 
Using (264) we find from (206) and (271) that the relative pressure perturbation is given by 

_ 2 A "^H^^ + i^BT(o) 

Op - On 3 ^ , „ 

"■B-t (0) 



The quantity A(t, a;) can now be rewritten as 



. 2^bT'(o) 

^ = 9,3, ^ - '^")- 
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Eliminating A(t, a;) from (281) with the help of (283) we find for the heat transfer to a perturbation from its sur- 
roundings 

T^„ysfl = ifcBTjo, [SSp - 5Sn) . (284) 

A density perturbation loses some of its heat energy to its surroundings if T^gysfl^ < 0, i.e. if 

5n{t,x)> I5pit,x). (285) 

If this condition is fulfilled, then a density perturbation may grow provided that its internal gravity is strong enough. 
The relative perturbation in the matter temperature and its time-derivative follow from (205), (264), (275), (278) and 
(280). We get. 



St^-IA St = 2HA, (286) 



so that for (5„ > we have St < 0, and the matter temperature perturbation increases slowly, St ^ 0. Combining 
(282) and (286), we arrive at the simple relation 

5p^5T + Sn, (287) 

i.e. both temperature fluctuations and perturbations in the particle number density contribute to perturbations in 
the pressure. Since tohc^ 2> kBT^o) during the era after decoupHng, it follows from (277), (279) and (286) that we 
may assume that 

Srit, x) « <5T(imat, a;), Srit, x) « 0, (288) 

to a very good approximation. Apparently the heat exchange is so efficient that the gas remains at a nearly constant 
temperature during the linear phase of the gravitational collapse. With (287) this implies that 

Sp{t,x) « ST{trm.t,x)+Sr,{t,x). (289) 

Thus, for a growing density perturbation the internal pressure will be increased. 
Expression (286) can be approximated by 

5„(t, x) - 6,{t, X) EE Ait, x) « _3 fcBT(o)(t) ^^ ^ '^<5T(t,nat, x), (290) 

2 tohc^ 10 

where we have used that mnc^ ^ fceTlo)- Note that (290) is a solution of (278), since T(o)(t) oc a^^(t), (277). 
Finally, using (289) the heat transfer (284) can now be rewritten as 

T^,,{t)st,{t,x) « iA:BT(„,(<)[3,5T(W,a;) - 25^{t,x)\. (291) 

This concludes the discussion of equation (201b). 

We now consider (201a). After substituting (271), (273) and (276) into the coefficients (202) it is found that 

6i=3i/, 62 = -|«;e,„,--|_, 63 = -3^, (292) 

considering that for a flat universe — V^. In the derivation of 63 we have used that <C 1, as follows from (266) 
and (273). For the evolution equation for density perturbations, (201a), this results in the simple form 

/^;2y2 \ 3i;2v2 

5e^iH5, ~ ^— -I- |k£,o) 5, = -^—ST{t^,t,x), (293) 
\ c ci J oca 

where we have used the result (290). Using (248) and (249) the evolution equation for the amplitude 5^{t,q) can be 
rewritten as 



5e + 3i/(t,nat) ( T^) 4 + i?'(tmat) 



Mm 



_8 

t \ 3 h f t 



^mat / ^ V ^mat 



_8 

^^t2u ^„2 / i \ 3 



(tmat)Mm ; '5T(imat,g), 

(294) 
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where we have incorporated (267)-(269), (273) and (277). The constant fi^ is given by 

Ws(iinat) 



1 



a(imat) -ff(imat) C 

Using the dimensionless time variable 



5 fcBr(0)(^mat) 
3 TTiH 



"7" = ~ ; ^mat ^ ^ ^ ; 
t-mat 



it is found from (268a) that 



d" 



1 y d" 
cimat / dr" 



2-H"(^mat) 



c"dt" 

Using this expression, equation (294) can be rewritten in the form 



10 



2 /4 (i^ 



S, = - 



d" 
dr" 



15 r8/3 



n = 1,2,. 



(>T{tma.t,q), 



(295) 



(296) 



(297) 



(298) 



where we have also used (290). In equation (298), a prime denotes differentiation with respect to r. 

It is of convenience for the numerical integration of equation (298) to express the dimensionless time variable t, 
(296), in the cosmological redshift z{t). Using (232) and (269), we get 



ait) 



3 

1\ 2 



,a(imat)/ V ^(^) + 1 

The integration will then be halted if either the time variable t has reached the value r^nd for which z = 0, i.e. 



(299) 



[^(tmat) + 1 



3/2 



(300) 



or when \S^{t, q)| = 1 has been reached for t < r^nd- In Section XV, we solve numerically the inhomogeneous equation 
(298). First, however, we consider the homogeneous part of the evolution equation (298). 



3. General Solution of the Evolution Equation 



In this section we calculate the exact solution of equation (298). To that end, we replace the independent variable 
T in this equation by the new independent variable 



so that 



3 Pm ' 



-7/3 



4,2 _-8/3_ 



dr s'^"'- dx' dT2 9^""^ dx 9^"^^ dx^' 



Using the variable (301), equation (298) can be rewritten in a much more tractable form 

d'^S^ 2d4 10' 



da;2 X dx 



1 - ^ 
x"^ 



STitmiit, q)- 



The general solution of this equation is 



AiJ 7{x)+A2J_7{x) 
2 2 . 



where Ai and A2 are the constants of integration and J±^{x) are the Bessel functions of the first kind: 



JMx) 



{x'^ — 15x) cos X — {6x^ — 15) sin x 



J 7 (x) = \ — (x — 15x) sinx + {6x — 15) cos a; 

"2 V TT 



,-7/2 



n-7/2 



(301) 
(302) 

(303) 

(304) 

(305a) 
(305b) 
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with the asymptotic expressions for a; — >■ 



'2 V TT 105 ^2 V TT X''^ 

Transforming back from a; to r, we arrive at the general solution of equation (298) 

3 



2 15 



Bi{q)J 7 (2^,„T 



-l/3^ 



■B2(q)J_7(2AimT 



-1/3^ 



r-1/2 _ 



1 



5t2/3 



ST{tmat,q), 



(306) 



(307) 



where Bi{q) and -62(9) are arbitrary functions. The first two terms are the solution of the homogeneous equation 
and the last term is the particular solution. The solution of the homogeneous equation follows from (5T(imat,'Z) = 
and has the property 



6T{t,x)=0 <S4> Sn{t,x) = 5^{t,x) — Sp{t,x), 



(308) 



where we have used (280), (283) and (286). Thus, if 5^ = (5„ then pressure perturbations counteract the growth or 
decay of a density perturbation in such a way that Sn — Sp and temperature perturbations in the matter do not occur. 
This is always the case in the standard Newtonian perturbation theory given by equation (357), since in this theory 
one has Sn = Sg. In the literature about linear cosmological perturbations, equation (278) does not exist, and the 
standard second-order evolution equation (357) is homogeneous. Since m^c^ 3> k-QT^g^, one is forced to take (5„ = S^. 
It follows from (284) and (308) that, in this case, the heat transfer to a perturbation is given by 



-t(0)*(l) 



(309) 



This implies that a density perturbation with (5„ > loses some of its internal heat energy to its surroundings, so 
that it may grow. The heat transfer T^o-^sfl^ from the perturbation to its surroundings is very small. In contrast to 
the standard theory, our treatise allows 5„ 7^ S^, which may result in a larger heat loss and, hence, a faster growth 
rate. This will be studied in Section XV. 

The evolution of density perturbations can be studied by imposing initial conditions (5e(imat,9) and 5^{tmat,q) on 
the general solution (307). Since the resulting expression is far too complicated, we investigate the evolution of density 
perturbations by solving equation (298) numerically in Section XV. In this section we only consider the two limiting 
cases of large-scale and small-scale perturbations. 

In the large-scale limit |q| — >■ 0, A — ^ c» (i.e. larger than the horizon, see Appendix D), it is found that, transforming 
back from t to t, 



2(5e(imat) 



<5e(imat) 



(310) 



Thus, for large-scale perturbations, the initial value (5T(^mat,Q) does not play a role during the evolution: large-scale 
perturbations evolve only under the influence of gravity. The solution proportional to t^/^ is a standard result. Since 
5^ is gauge-invariant, the standard non-physical gauge mode proportional to t^^ is absent from our theory. Instead, 
a physical mode proportional to t~^/'^ is found. This mode has also been found by Bardeen [12], Table I, and by 
Mukhanov et al. [13], expression (5.33). In order to arrive at the t~^/^ mode, Bardeen has to use the 'uniform Hubble 
constant gauge.' In our treatise, however, the Hubble function is automatically uniform, without any additional gauge 
condition, see (25). 

In the small-scale limit A — > or, equivalently, \q\ -> 00, we find, transforming back from r to t, 



Seit,q) ~ -Uritmat^q) 



t 

^mat 



-1/3 



<5T(imat,g) +<5e(imat,q) 



cos 



2/^ni - 2/i 

n 



-1/3- 



(311) 



Thus, small-scale density perturbations oscillate with a decaying amplitude which is smaller than unity so that the 
non-linear regime will never be reached. 

For scales between the two extremes discussed above the growth of the perturbations can be considerable as we will 
show in Section XV. 
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XIV. STAR FORMATION: BASIC EQUATIONS 



This section is a preparation for the numerical solution of the evolution equation (298). We express the constant 
/j,„i and the star mass Af(imat) hi the observable quantities z(tdec), 'H{tp) and T(o)(tdec), the initial redshift z(tmat)) 
where imat is some initial time (270) at which a density perturbation starts to contract, and the physical dimensions 
Amat of a density perturbation. Finally, we study the influence of the particle mass m and the initial time fmat on the 
mass of a star. 



A. Initial Values and the Mass of a Star 



Writing q — 2tt/X, where Aa(iniat) is the physical scale of a perturbation at time ii„at and Aa(ip) = A is the physical 
scale as measured at the present time tp, we get for (295) 



^ /5 fcBr(0)(^mat) , - \ Ci^O\ 

A^m = T 7777 ^\ Q 2 ' -^^at = Aa(t,natj- (312) 

Amat H{t^^t) V 3 miiC^ 



The Hubble function i?(imat) follows from (232), (268a) and (269). We find 

i?(imat) = H{tp)[z{t^,t) + 1]^/'. (313) 

From (277), we get 



^(imat) + 1 



2 



T,,(.mat)=r<,(W)^^g^J . (314) 

Using (313) and (314), expression (312) can be rewritten as 



Amat H{tp) [^(tdcc) + 1] V^(^mat) + 1 V ^ "^H ' 

where we have used (54). With (315) we have expressed /im in the observable quantities ■H(ip), 2;(tdoc) and T(o)(tdcc)- 
The matter temperature just after decoupling, T(o)(tdcc), is given by (265). Using also (231) we get 

512.0 . 

Mm = , =, Amat m pC, (316) 

Amat V ^(^mat) + 1 

where we have used that 1 pc = 3.0857 x 10^^ m (1 pc = 3.2616 ly). 

The mass Af (imat) of a spherical density perturbation with radius |Amat is given by 

47r 3 

A^(tmat) = (5 Amat) ^^(0) (imat)mH, (317) 

where toh is the proton mass. The particle number density can be calculated from the value at the end of the radiation- 
dominated era. By definition, at the end of the radiation-domination era the matter energy density n(o)TOHC^ equals 
the energy density of the radiation. Hence 

?^(0)(icq)?7T-HC^ = aBT'(„)^(tcq). (318) 

Since n(o) oc and T^a■^~^ cx a~^, we find, using (232), the particle number density at time imat 

03^(0)^ (ip) 



^(0)(tmat) = W-^) + 1] [^(^-at) + 1] . (319) 



Combining (317) and (319), we get 



M(imat) = ^ (iAmat)' ^^Z^^^ + l] + l]^. (32O) 

O C 
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Using that one solar mass is 1.98892 x 10^'' kg, we find from (231) 

M{U,t) = 1.141 X 10-«ALt[z(t„,at) + 1]'M0, A,„at in pc (321) 

This expression will be used to convert the scale Amat of a perturbation, which starts to contract at a redshift of 
z{tniat), into its mass. The latter is expressed in units of the solar mass. It should be stressed here that the numeric 
factors in expressions (316) and (321) hold true only for a fluid consisting of protons. In Section XIV B we investigate 
the influence of the particle mass on the mass of a star. 

In order to solve equation (298), we need the initial values Seitmat, q), S'^{tmat,(l) and dritmat, q)- From (290) it 
follows that, since \5T{t,q)\ < 1, we must have 

Mt,q)-S,it,q)\ « 2|^^(t,„,„q)|MkW < g) | Ml^ll^ ^ 4.10 X lO^i" |<5T(imat, q)| , (322) 

implying that 

5.a{t,q)^5e{t,q), t>tdec- (323) 

We take, however, 5„(tmat:'7) 7^ '^e(^mat5 9)- The case (5„(tmat:'7) = '^e(^mat5 9) has been discussed in Section XIII C 3 
on the general solution of the evolution equation (298). 

An initial value for the relative matter temperature perturbation <^T(^mat,Q) can be found as follows. A weak 
condition for growth is that a perturbation must lose some of its internal heat energy, i.e. T^gysfl^ < or, equivalently, 
(285). A larger growth can be achieved if we take the initial values such that the heat loss of a density perturbation 
is larger than given by (309), i.e. 

- kBT^„,{t)Sn{t,q) > ifcBT(o,(i)[3<5T(imat,g) ^ 2(5„(i,q)] ^ SriUauq) < 0, (324) 

where we have also used the combined First and Second Laws of thermodynamics in the form (291). The criterion (324), 
yields a positive source term of the evolution equation (298) and, therefore, a larger growth than the homogeneous 
equation. 

Finally, there are no observations of the growth rates of density perturbations. Therefore, we take the initial growth 
rates equal to zero, i.e. 

^e(^mat, q) = <5ji (tmat , Q) = 0, ^doc < ^mat < • (325) 

In the next two subsections we investigate the influence of both the particle mass and the initial time of star 
formation on the mass of a particular star. In Section XV we present the results of our new perturbation theory for 
structure formation. 



B. The Influence of the Particle Mass on the Mass of a Star 



Up till now we have assumed that the matter content of the universe after decoupling consists of protons. However, 
WMAP observations suggest that there may exist a considerable amount of, as yet unknown, dark matter (dm) particles. 
In this section we investigate the influence of the mean particle mass on the evolution of density perturbations. We 
assume that after decoupling, a gas consisting of particles with mean mass ifi, given by {a{t) > 0) 

m = a(imat)'7iH, (326) 

can be described by an equation of state of the form (264) with the mass mn replaced by (326). As a crude estimate 
for m, we take the mean mass of the baryons and the dm particles, i.e. 

Q^(^) = — f m ' (327) 

where Qhar and 57dm are the baryon and dark matter particle densities respectively, given in units of the critical 
density, see the Friedmann equation (155). 

By /im, we denote the parameter fi^ (315) in which mjj is replaced by o:(t^^t)mYi. Density perturbations in a gas 
with mean particle mass m evolve in exactly the same way as perturbations in a gas of which the particles have mass 
77iH if /im = jj-ni- in this case we have 6e{t, q) = S^{t, q). From (315) it follows that 
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where Amat is the scale of a density perturbation in a gas with mean particle mass fh. Using (317), (326) and (328) 
we find 

M(imat) = ^p^, (329) 
V°^(*mat) 

where -/Vf(ii„at) refers to a perturbation mass in a gas with mean particle mass rh. In the derivation of (329) we have 
assumed that n(o)(iinat) = 'T^(n)(imat), i-e. the particle number density is for a mean particle mass rh the same as for 
a particle mass toh- We thus have found that perturbations in a gas consisting of particles which are heavier (i.e. 
Q^(^inat) > 1) than protons the collapse takes place at a smaller total mass than perturbations in ordinary matter. In 
other words, heavier particles yield lighter stars. We will illustrate (329) with two extreme cases, namely a mixture 
of heavy wimps (cdm) with a mass of mcdm ~ 70 mn [38] and baryons, and a mixture of ordinary matter and 
hot dark matter (hdm) with a mass of mhdm = 1.5eV/c^, as suggested by Nieuwenhuizen et al. [5, 40]. Using that 
Toh = 0.938 GeV/c^, we can calculate with the help of (157) and (327) the mean particle mass and, hence, a. Assuming 
that ribar(i)/f^dm(i) ~ f^bar(ip)/i^dm(ip) for both typcs of particlcs and for t > t^ec, we find a{t) ~ a{tp), so that 

1/2 1/2 

•^cdm ^ 0.13 and Oj^^j^ ~ 2.4. This implies that if the dark matter particles are heavy WiMPs then the stars formed 
are much lighter than the stars formed in a universe filled with baryons only. On the other hand, if light neutrinos 
are the dark matter (hdm), then the stars formed will be heavier. 



C. The Influence of Initial Time on the Mass of a Star 



In this section we show that density fluctuations which start to contract at late times yield stars with smaller masses 
than early density perturbations. To show this, we consider the evolution equation (298). 

It follows from equation (298) that perturbations starting to grow at t,„at > ^dec and obeying the same initial 
conditions as fluctuations starting to contract at idcc; i-e. 

6e{tynat,x) ^ Se{tdcc,x), S'^{t^s^t , x) ^ 5'^{tdcc, x) , Srit^at, x) ^ Sritdcc, x) , (330) 

evolve in exactly the same way, provided that 

A*m(^mat) = A*m(idcc), (331) 

where we have written /im(io) to denote the value of the parameter /im taken at time t ~ tQ. Using the relation 
(331) we can relate the masses of density perturbations starting at different times. Replacing the initial time i,„at in 
equation (298) by the initial time idecj we get for the parameter (315), at t = t^ec 



„ f, N 27r 1 /5 fcBT(o)(^doc) .„„„^, 

Atm (Idee) = T -JJ^ \ / o ■ (^32) 

^dccH(ip)[z(tdcc) + l] ^ V 3 "^H 



Equating (315) and (332), we arrive at 



Amat — Adoc-l/ ! , -, j (333) 

y ^(imat) + 1 

i.e. the evolution of a density perturbation starting at t ~ t^cc with scale Adcc is exactly equal to the evolution of a 
density perturbation starting at i = imat with scale Amat, if and only if the relation between the scales obeys (333). 
Since n(o) oc cx (z + l)'^, we have 

n^o){tmat) = »(o)(W) [ "I" ) . (334) 



zitdcc) + 1 



Using these expressions and (317), we flnally arrive at 



M(W) = A/(idec) (^Tr4^) ' ■ (335) 

V 2(^dec) + 1 / 

This expression relates the masses of density perturbations, which start to evolve at different times idee and tmat with 
the same initial values (330) and (331). Consequently, stars formed from late time fluctuations have smaller masses. 
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XV. STAR FORMATION: RESULTS 



The standard cosmological theory of small perturbations is characterized by 6n — S^, (308). This is, however, too 
restrictive. Although \Sn ~Se\ is very small, it need not be zero, as follows from (322). Since mac^ ^ fceTJo), it follows 
from (282) and (286) that the quantities St, Sp and 5„ are not, beforehand, confined to small values: small changes 
in A = (5„ — (5e may lead to large changes in St, Sp and Sn- These quantities must only fulfill the linearity conditions 
\St\ < 1, \Sp\ < 1 and \Sn\ < 1. The fact that Sn need not be exactly equal to S^ is paramount for the formation of 
structure in the universe. 

In this section we will solve the evolution equation (298) numerically. To that end we use the differential equation 
solver Isodar with root finding capabilities, included in the package deSolve, which, in turn, is included in R, a system 
for statistical computation and graphics [41]. We investigate star formation which starts at cosmological redshifts 
z = 1091 and z ^ 1. 



A. Star Formation starting at z — 1091: Population III Stars 



At the moment of decoupling of matter and radiation photons could not ionize matter any more and the two 
constituents fell out of thermal equilibrium. As a consequence, the pressure drops from a very high radiation pressure 
P = ^clbT^ just before decoupling to a very low gas pressure p = nk^T after decoupling. This fast transition from a 
high pressure epoch to a very low pressure era may result in large relative pressure perturbations. In this subsection 
we study the influence of these relative pressure perturbations on the formation of stars. 



1. Initial Values 



In order to integrate equation (298), we need initial values. The initial value S^(t^cc,<i) is related to the relative 
perturbation (5T^(idcc,Q) in the background radiation. Using (205) and (261), we get at the end of the plasma era 



aBT(o)^(tdec) 



[5e {tdec ,q)-Sn (<dec , Q)] + Sg (tdec , 



Using that n^o) oc a ^ and TJoj^ cx a ^, we get from (232) and (318) 



'T-(0)(idec)wHC^ = aBr"™ (tdoc)- ^ 



z{t 



dec 7 



1' 



so that (336) can be rewritten as 



5T^(tdec,q) = ^ 



zit. 



cqj 



-y [Se{tdec, q) - Snitdec, ?)] + '^^(idec, Q) 



_z{tdcc) + 

In order to eliminate 5„(tdcc, q), we use that at the end of the plasma era we have, according to (263), 

^e{tdec,q) 



where, using (261) and (337), 



w{td 



Sn{tdec, q) - 



P(0) (^dcc) 



1 + w(idec) 
z{tdcc) + 1 



0, 



£(0)(tdec) 3z(tcq) + 3z(tdcc) + 6 



0.085. 



Combining (339) and (340), we find 



. ,, ^ J,, , 3z(tcq) + 3z(tdcc) + 6 
iz{tcc^) + 4z(tdcc) + 7 



We can now rewrite (338) as 

ST.yitdec,q) = Ss{tdec,q) 



z(^eq) + -g(^dcc) + 2 
30(teq) +4z(tdcc) +7 



0.92 (5, (tdec, q). 



0.31 Ss{tdec,q) 



(336) 



(337) 



(338) 



(339) 



(340) 



(341) 



(342) 
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Star Formation starting at z = 1091 




I I I I I I I I I I I L 

le+01 le+02 le+03 le+04 le+05 le+06 

Star Mass in Units of the Solar Mass 



Figure 1: The curves give the redshift at which a linear perturbation in the particle number density with initial values 
Sn{tdec,Q) ~ 10~^ and 5^tdec,q) = starting to grow at an initial redshift of z{tdcc) = 1091 becomes non-linear, i.e. 5„ ~ 1. 
During the evolution we have 5p{t,q) — ST{tdec,q) + 5„{t,q). The numbers at each of the curves are the initial relative 
perturbations in the matter temperature Sritdec, q)- For each curve, the maximum is at 3.4 x 10^ Mq. If cdm is present then 
the peak mass is at 4.5 x 10^ Mq. 

From the WMAP observation (231f) we find, using (342), 

Seitdcc, q) ~ 3.3 X 10-^ j;(idcc, q) ~ 0, (343) 

where we have assumed, for lack of observations, that during the decoupling of matter and radiation the growth rate 
5'^ is very small, see (325). With (343) the condition (323) implies that 

5„(idcc, q) ~ 3.3 X 10-^ S:,{Uec, q) ~ 0, (344) 

i.e. WMAP observations demand that, just after decoupling, also the relative particle number density perturbations 
are very small. 



As we have shown, see (288), the relative matter temperature perturbation St is very nearly constant for a contract- 
ing or expanding density perturbation. From the fact that the pressure is given hy p = nk-eT and (344) it follows that 
large pressure perturbations can, just after decoupling, only be realized by large matter temperature perturbations, 
hi our calculations we take, therefore, |(5T(^dcc: 3;)| in the range 0.5%-lG%. With (287) and (344) this implies that, 
initially. 



5p{tdcc, q) ~ Sritdec, q)- 



(345) 
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In other words, the perturbation in the pressure is, initially, mainly determined by a perturbation in the matter 
temperature. In view of (324) and (344) we have to take for the initial value of the relative matter temperature 
perturbation 

ST{tdcc,q) < 0, (346) 

in order to find growing density perturbations. We now have gathered the necessary ingredients to integrate the 
evolution equation (298) numerically. 

2. Results 

Figure 1 has been constructed [41] as follows. For each choice of SxitdecQ) we integrate equation (298) for a large 
number of values for Adec, using the initial values (343) and (344). The integration starts at t = i/tdec = 1, i-e. at 
z = z(tdoc) and will be halted if either z = [i.e. r = Tend, (300)[ or 5e{t,q) = 1 has been reached. One integration 
run yields one point on the curve for a particular choice of Adec if ^eitjQ) — 1 has been reached for z > 0. If the 
integration halts at z = and S^{t, q) < 1, then the perturbation belonging to that particular Adec has not yet reached 
its non-linear phase today, i.e. at t = 13.7Gyr. On the other hand, if the integration is stopped at Se{t,q) = 1, then 
the perturbation has become non-linear within 13.7 Gyr. In Figure 1 we have used the star mass M{tdec), expressed 
in solar masses, instead of the scale Adec of a density perturbation. To that end we have used expression (321). 

The above described procedure is repeated for 6T{tdcc, q) in the range —0.005, —0.01, —0.02, . . . , —0.1. During the 
evolution, the relative pressure perturbation evolves as 

Sp{t, q) = Sritdec, q) + Sn{t, q). (347) 

The fastest growth is seen for perturbations with a mass of approximately 3.4 x 10'^ M©. This value is nearly 
independent of the initial value of the matter temperature perturbation ^^(idccj q)- Even density perturbations with 
an initial matter temperature perturbation as small as SxitdccTq) = —0.5% reach their non-linear phase at z = 0.13 
(T(o)^ = 3.1 K, t — 11.5 Gyr) provided that its mass is around 3.4 x lO'^M©. Perturbations with masses smaller than 
3.4 X 10"^ M0 reach their non-linear phase at a later time, because their internal gravity is weaker. On the other hand, 
perturbations with masses larger than 3.4 x 10'^ Mq have larger scales so that they cool down slower, resulting also 
in a smaller growth rate. Since the growth rate decreases rapidly for perturbations with masses below 3.4 x 10'^ Mq, 
the latter may be considered as the relativistic counterpart of the Jeans mass. 

Figure 1 has been calculated for a baryonic cosmic fluid, without CDM or hdm. If CDM is present, then it follows 
from (329) that the graphs shift towards smaller masses such that the Jeans masses are at 4.5 x 10^ M©. For hdm 
the peaks are at 8.2 x 10^ M©. 

Finally, using (321), we find that all star masses in Figure 1 start to contract at decoupling from density per- 
turbations with diameters less than 52 pc, which is much smaller than the particle horizon size (Appendix D) 
dnitdcc) = 349 kpc. 

B. Star Formation starting at z; = 1 

At z = 1 the interstellar gas has been diluted so much that star formation is not possible anymore in the interstellar 
gas with only the negative matter temperature perturbation STitmut, q) as the driving force. Consequently, late time 
star formation can only take place in regions which have a higher density with respect to the intergalactic space. 
Therefore, late time star formation takes place mainly within galaxies. As can be seen from Figures 2 and 3, the 
initial density perturbation should be at least of the order of 5„ « (5e « 0.70, in order to yield eventually a gravitational 
collapse for z > 0. 

We have considered two extreme cases of star formation starting at z(tniat) — 1 or, equivalently, tmat ~ 4.8 Gyr. 
In both cases, we assume that the initial growth rate vanishes, (325). In the first case, depicted in Figure 2, we 
have chosen Srit^at, q) = —1, the smallest value that a relative perturbation can have without violating the linearity 
conditions. During the evolution of a density perturbation, we have, according to (289), 

Sp{t,q) = -l + 5n{t,q), (348) 

so that during the evolution the relative pressure perturbation is negative. This is the most favorable situation for a 
perturbation to grow: pressure cooperates with gravitation. The most conspicuous feature in Figure 2 is the sharp 
lower limit of star formation. Density perturbations with masses below 0.2 M©, do not become non-linear before 
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Figure 2: Star formation for &T{t, q) = — 1- The curves give the redshift at which a linear density perturbation starting to grow 
at an initial redshift of z(tmat) ~ 1 becomes non-linear, i.e. 5„ ~ (5e ~ 1. The numbers at each of the curves are the initial 
relative density perturbations (5n(imat,q) ~ Ss{tnia.t,q)- During the evolution we have Sp{t,q) = —l + 5„{t,q). For masses in 
excess of 100 M©, the growth becomes independent of the scale of a perturbation, i.e. the growth is proportional to t^^^. 



z = is reached. Apparently, the gravitational field is. for masses below 0.2 M©, too weak to collapse in due time 
and, eventually, become a star. Another notable characteristic of star formation is that it is slow: even if the initial 
density perturbation is as high as (5„(i,„at,Q) — 0.999 it takes, for a 100 M0 perturbation, approximately 214Myr to 
reach the value 6n{t,q) = 1. 

The second case, characterized by 5T(^mat,'7) = and summarized in Figure 3, is the most detrimental for star 
formation, because during the evolution of a density perturbation we have 

Spit,q)^5nit,q), (349) 

so that the pressure opposes the contraction of a density perturbation. However, if the internal gravity of a per- 
turbation is strong enough, then gravity will overcome the pressure and the perturbation will eventually collapse to 
form a star. Due to the counteracting pressure, the lower limit of star formation is much larger than in the case of 
cooperating pressure, namely 0.8 M©. 

Both extreme cases have one common characteristic. For perturbations with masses larger than 100 , the collapse 
time is nearly the same, as can be seen from Figures 2 and 3. Apparently, the internal gravitational field is so strong 
that opposing or cooperating pressure perturbations do not play a role anymore. 
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Figure 3: Star formation for Srit, q) = 0. The curves give the redshift at which a linear density perturbation starting to grow 
at an initial redshift of z(tmat) = 1 becomes non-linear, i.e. 5n = 5^ ~ 1. The numbers at each of the curves are the initial 
relative perturbations 5n(tmat, q) = ^^(tmat, q)- During the evolution we have 5p{t, q) = 5n{t, q) — 3^(1, q). For masses in excess 
of 100 Mq, the growth becomes independent of the scale of a perturbation, i.e. the growth is proportional to t^^^ . 



XVI. STANDARD NEWTONIAN THEORY OF COSMOLOGICAL PERTURBATIONS 

The new evolution equations (245a) and (293) are different from their standard counterparts (350) and (357) 
respectively. In this section we explain the differences and we show why equations (350) and (357) should not be used 
anymore in the study of cosmological density perturbations. 

Padmanabhan had already in 1993 the supposition that the Newtonian theory of cosmological density perturbations 
is questionable. On page 136 of his textbook [39] he states: 

To avoid any misunderstanding, we emphasize the following fact: it is not possible to study cosmology 
using Newtonian gravity without invoking procedures which are mathematically dubious. 

Padmanabhan does not conclude, however, that the Newtonian theory of cosmological perturbations is incorrect, since 
he subsequently states: 

However, if we are only interested in regions much smaller than the characteristic length scale set by the 
curvature of space-time, then one can introduce a valid approximation to general relativity. 

Up till now, this point of view is considered as standard knowledge. This is due to the fact that the standard equation 
(357) of the Newtonian theory for small-scale perturbations is, in the low velocity limit Ws/c — > 0, similar to the 
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relativistic equation (362) derived from the General Theory of Relativity. This similarity has led to much confusion 
in the literature (see, for example, [32-34, 42, 43]), as we will now explain. First, we remark that gauge problems, 
although time and space coordinates may be chosen freely according to (220), do not occur in the Newtonian theory 
of gravity. This is because in the Newtonian theory itself the universe is static, i.e. £(0) = and n(o) — 0, implying 
with (219) that both e^r-, and tiji, are independent of the choice of a system of reference. Now, the reasoning in the 
literature is as follows. The quantity 6 = ea^/Sfo) occurring in the Newtonian equation (357) is, according to the 
standard knowledge, independent of the gauge choice, since this equation is derived from the Newtonian theory, in 
which gauge problems concerning perturbed scalar quantities do not occur. Furthermore, the Newtonian equation 
(357) is valid only for density perturbations with scales smaller than the horizon size. By virtue of the resemblance 
of the Newtonian equation (357) and the relativistic equation (362) it is, therefore, put forward that 

a gauge dependent quantity, such as 6 = e(i)/e(o), which is initially larger than the horizon, becomes 
automatically gauge-invariant as soon as the perturbation becomes smaller than the horizon. 

This viewpoint is, however, incorrect as we will now show in detail. Firstly, we remark that the universe cannot be 
static in the non-relativistic limit (see Section XII) of the General Theory of Relativity, since — >■ violates the 
Einstein equations (214). Secondly, it has been shown in Section XII that there remains some gauge freedom in the 
non-relativistic limit, namely the freedom to shift time coordinates, x*^ — >■ — ip, and the freedom to choose spatial 
coordinates, a;* — >■ a;* — x*(^)- This coordinate freedom is a well-known and natural property of Newtonian physics 
and, therefore, should follow from a relativistic perturbation theory in the low velocity limit Vg/c — > 0. Thirdly, 
we will show that, just because of the resemblance of the relativistic and Newtonian equation, the gauge function 
tp also occurs in the solution (359) of the Newtonian equation (357). Consequently, the standard equation of the 
Newtonian theory has no physical significance: due to the appearance of the gauge function -0 in one of the two 
independent solutions, one cannot impose initial conditions to arrive at a physical solution. From the fact that, in 
the non-relativistic limit of the General Theory of Relativity, the universe is not static, combined with the freedom 
in time coordinates implies that if a quantity, such as d = e(^/e^o■n is gauge dependent in the General Theory of 
Relativity it is also gauge dependent in the non-relativistic limit, as follows from (219). The important conclusion is, 
therefore, that a quantity which has no physical significance outside the horizon, does not become a physical quantity 
inside the horizon. This conclusion is consistent with the facts that, in the non-relativistic limit of our perturbation 
theory, the gauge-invariant quantities ef'), (225), and nfly, (228), survive, while the gauge dependent quantities £(1) 
and 71(1) disappear from the scene, as we have shown in detail in Section XII. Hence, with efl^ and nfl^ surviving in 
the non-relativistic limit, there is indeed no gauge problem in the Newtonian theory of gravity. 

The fact that linear perturbation theory is plagued by the gauge solutions (356) and (363) has already been pointed 
out by Lifschitz [6, 22, 23] as early as 1946. In 1980, thirty-four years later. Press and Vishniac ]11] called attention 
to the same issue. In spite of these warnings, the standard equations (350) and (357) are still ubiquitous in the 
cosmological literature. Apparently, the cosmological gauge problem is quite persevering. 

In the next two subsections, we will elucidate the gauge problem. We go one step further than Lifschitz and Press 
and Vishniac. These researchers have shown that only for large scales the solutions of the standard equations are 
infected by spurious gauge modes. We show, in addition, that the solutions (352) and (359) of the standard equations 
(350) and (357) contain the gauge function ip{x), independent of the scale of a perturbation. Consequently, the 
important conclusion must be that 

the Newtonian theory of gravity is not suitable to study the evolution of cosmological density perturbations. 

In order to show that the Newtonian theory of cosmological perturbations is invalid, we consider a flat {k — 0) flrw 
universe with a vanishing cosmological constant (A ~ 0) in the radiation-dominated era and the era after decoupling 
of matter and radiation. 



A. Radiation-dominated Universe 

The standard equation for the density contrast function 5 which can be found, for example, in the textbook of 
Peacock [44], equation (15.25), is given by 

^ + 2i/^- Q^ + |Ke(o)^ ^ = 0. (350) 

Compare this equation with equation (245a) for the gauge-invariant contrast function d^- Equation (350) is derived 
by using special relativistic fluid mechanics and the Newtonian theory of gravity with a relativistic source term. In 
agreement with the text under equation (15.25) of this textbook, the term —^W^5/a^ has been added. The same 
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result, equation (350), can be found in Weinberg's classic [35], equation (15.10.57) with p = and Vs = l/\/3- Note, 
that equation (350) cannot be derived from the General Theory of Relativity. 

Using (240), (241), (248), (251) and (252), we can rewrite equation (350) in the form 



^" + V+(g-^)^ = 0, (351) 



where the constant /ij. is given by (254). A prime denotes differentiation with respect to r, (251). The general solution 
of this equation is found to be 

d{T,q) = ^^^J2{^^rV^)+i^iq)n^l^.H{t,,^)Y4^^,V^), (352) 

where Ci{q) and ip{q) are arbitrary functions (the integration 'constants') and J,j(x) and Y^{x) are Bessel functions 
of the first and second kind respectively. The factors S/fi^ and 7r/i^i/(trad) have been inserted for convenience. Thus, 
the standard equation (350) yields oscillating density perturbations with a decaying amplitude. 

For large-scale perturbations (|q| — > or, equivalently, /^r ^ 0), the asymptotic expressions for the Bessel functions 
J2 and Y2 are given by 

J2(/irV^) « 'fr, Y2{fi,V^)^-^T-\ (353) 

^ ' 8 TT/iJ 

Substituting these expressions into (352) , it is found for large-scale perturbations that 

S{t) = CiT - AH{t,^A)^T-\ (354) 

where we have used that Ci(|q| ^ 0) = Ci and '0(|q| ^ 0) = become constants in the large-scale limit. Large-scale 
perturbations can also be obtained from the standard equation (350) by substituting = 0, i.e. 

(5 + 2i7(5- |k£(o)(5 = 0. (355) 

The general solution of this equation is, using (351) with /if = O, given by (354). Thus far, the functions Ci{q) and 
il>{q) are the integration 'constants' which can be determined by the initial values (5 (trad i?) ^^nd (5(trad7 9)- 

However, equation (355) can, in contrast to equation (350), also be derived from the General Theory of Relativity: 
see the derivation in Appendix E. As a consequence, equation (355) is found to be also a relativistic equation, implying 
that the quantity d — e^^/e^o) is gauge dependent. Therefore, the second term in the solution (354) is not a physical 
mode, but equal to the gauge mode 

'5gauge(^) = = -4i/(irad)V'T-\ (356) 

£(0) 

as follows from (2a), (239b) and (240a). Consequently, the constant ip in (354) and, hence, the function ■ip{q) in 
(352) should not be interpreted as an integration constant, but as a gauge function, which cannot be determined 
by imposing initial value conditions, see Appendix C for a detailed explanation. Thus, the general solution (352) of 
the standard equation (350) depends on the gauge function ip{q) and has, as a consequence, no physical significance. 
This, in turn, implies that the standard equation (350) does not describe the evolution of density perturbations. 

Here the negative effect of the gauge function is clearly seen: as yet it was commonly accepted that small-scale 
perturbations in the radiation-dominated era of a flat flrw universe oscillate with a decaying amplitude, according 
to (352). The treatise presented in this article reveals, however, that small-scale density perturbations oscillate with 
an increasing amplitude, according to (260). This is the real behavior of a small-scale density perturbation. 

B. Era after Decoupling of Matter and Radiation 

The standard perturbation equation of the Newtonian theory of gravity is derived from approximate, non-relativistic 
equations. It reads 

d + 2HS- i-^-^ + lKs^o^S = 0, (357) 
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Star Formation starting at z = 1 : Newtonian Theory 
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Figure 4: Star formation according to the standard equation (357) with 5T{t,q) = 0. The curves give the redshift at which a 
linear density perturbation starting to grow at an initial redshift of zitmat) ~ 1 becomes non-linear, i.e. 5n = = Sp ^ 1. The 
numbers at each of the curves are the initial relative density perturbations 5n(tmat,q) ~ Se{tina.t,q) ~ 5p(tmat,q). During the 
evolution we have 5p{t,q) — Sn{t,q) = Ss{t,q). For masses in excess of 100 Mq, the growth becomes independent of the scale 
of a perturbation, i.e. the growth is proportional to t'^^^. 



where is the speed of sound. (See, for example, Weinberg [35], Section 15.9, or Peacock [44], Section 15.2.) Compare 
this equation with equation (293) for the gauge-invariant contrast function S^- Just as equation (350), the standard 
equation (357) cannot be derived from Einstein's gravitation theory. 

Using (248), (268), (269), (296) and (297), equation (357) can rewritten in the form 



S" 



3t 



-S' 



4 f^l 



9 r8/3 3t2 



0, 



(358) 



where the constant /im is given by (295). A prime denotes differentiation with respect to t, (296). Tlie general solution 
of equation (358) is found to be 



^l(9)y^J_5(2AimT-l/3)-^7A(q)i?(t,„at),/^J 5(2/i„,T-l/3) 



-1/6 



(359) 



where Di{q) and ijj{q) are arbitrary functions (the 'constants' of integration) and J^jj[x) is the Bessel function of the 
first kind. The factors |-\/7r/i^ and ^i?(i,„at)-\/ '''/Mm h^^ve been inserted for convenience. 



57 



We now consider large-scale perturbations characterized by = (i.e. |<7| 0) or perturbations of all scales 
in the limit Vs/c — > 0. Both limits imply /Xm — > 0, as follows from (295). The asymptotic expressions for the Bessel 
functions in the limit fi^i ^ are given by 




(360) 



Substituting these expressions into the general solution (359), results in 

5{t) = D,T^^^ - 3ff(i„at)V'r-i. (361) 

where we have used that i'ldq] 0) = Di and V'd?! — >■ 0) = '0 become constants in either the large-scale limit 
|q| or in the limit v^/c — > 0. In the limit /i^ — >■ 0, equation (357) reduces to 

S + 2HS-^Ke(„^6^ 0. (362) 

Using (358) with = 0, we find that the general solution of equation (362) is given by (361). Thus far, the functions 
Di{q) and 'ip{q) are the integration 'constants' which can be determined by the initial values (5(<inat, o) and (5(tinat, o)- 
However, equation (362) can, unlike equation (357), also be derived from the General Theory of Relativity, and is, 
as a consequence, a relativistic equation: see Appendix E for a derivation. In this case, however, it is based on the 
gauge dependent quantity S = e^l■)/e^o)■ As a consequence, the second term of (361) is equal to the gauge mode 

Wc(t) = = -3i/(t„,at)^r-\ (363) 

£(0) 

as follows from (2a), (267b) and (268a). Therefore, the constant tp in (363) and, hence, the function tp{q) in (359), 
should not be interpreted as an integration constant, but as a gauge function, which cannot be determined by imposing 
initial value conditions (see Appendix C). Since the solution (359) of equation (357) depends on the gauge function 
ip{q) it has no physical significance. Consequently, the standard equation (357) does not describe the evolution of 
density perturbations. 

Again, we encounter the negative effect of the gauge function: up till now it was commonly accepted that for small- 
scale density perturbations (i.e. density perturbations with wave lengths much smaller than the particle horizon. 
Appendix D) the Newtonian theory suffices and gauge ambiguities do not occur and that the evolution of density 
perturbations in the Newtonian regime is described by the standard equation (357). The treatise presented in this 
article reveals, however, that in a fluid with an equation of state (264), the evolution of density perturbations is 
described by the relativistic equation (293), for small-scale as well as large-scale perturbations. In fact, equation (293) 
explains the formation of stars in the universe. 

For pedagogical reasons only, we have calculated the 'formation of stars,' starting at z = 1, with the help of the 
standard equation (357) for which the conditions (308) hold true. The result is depicted in Figure 4. This figure 
should be compared with Figure 3, which is calculated with the help of the new equation (293) under the conditions 
(308). The standard Newtonian theory underestimates the internal gravitational field of a perturbation by a factor 
of I with respect to our perturbation theory. Compare the gravitational fields ^k£(o) in equation (357) and its 
relativistic counterpart (362) with |k£(o) in the new equation (293). As a result, the standard theory predicts a lower 
limit of 1.7M0 for star formation, whereas our theory yields a lower limit of O.SMq under the conditions (308). As 
a consequence, our Sun could not exist at all according to the standard theory. This fact can be considered as an 
'experimental proof that the standard Newtonian theory has indeed no physical significance. 



Appendix A: Equations of State for the Energy Density and Pressure 

We have used an equation of state for the pressure of the form p = p{n,e). In general, however, this equation of 
state is given in the form of two equations for the energy density e and the pressure p which contain also the absolute 
temperature T: 

e = e{n,T), p^p{n,T). (Al) 

In principle it is possible to eliminate T from the two equations (Al) to get p — p{n, e), so that our choice of the form 
p — p{n,e) is justified. In practice, however, it may in general be difficult to eliminate the temperature T from the 
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equations (Al). However, this is not necessary, since the partial derivatives and p„ (78), the only quantities that 
are actually needed, can be found in an alternative way. From equations (Al) it follows 

From (A2b) it follows that the partial derivatives (78) are 

..^(i)m. (A3b, 



(A4a) 

Upon substituting the expressions (A4) into (A3), we find for the partial derivatives defined by (78) 

- (II ^ (SI- (11(11 



From (A2a) it follows 



' ^ -{^] ■ (A4b) 



(A5b) 



where e and p are given by (Al). In order to calculate the second-order derivative p„„ replace p in (A5a) by p„. For 
Pee replace p in (A5b) by p^. Finally, for pen = Pna replace p in (A5a) by p^ or, equivalently, replace p in (A5b) by p„. 

Appendix B: Derivation of the Manifestly Gauge- invariant Perturbation Equations 

In this appendix we derive the perturbation equations (197) and the evolution equations (201). 



1. Derivation of the Evolution Equation for the Entropy 

With the help of equations (195a) and (195b) and equations (153)-(154) one may verify that 

In view of (174) one may replace £(i) and n(i) by efl^ and nfly Using (177) yields equation (197b) of the main text. 

2. Derivation of the Evolution Equation for the Energy Density Perturbation 

We will now derive equation (197a). To that end, we rewrite the system (195) and expression (196a) in the form 



£(1) - 


t- aiiEd) - 


1- Q;i2n(i) - 


1- ai3i9(i) 


-t-ai4^i?(i)|| 


= 0, 


(B2a) 


fiii) - 


1- a2ie(i) - 


f- Q;22n(i) - 


f- ^23^9(1) 


+ a24^i?(i)|| 


= 0, 


(B2b) 




1- asie,!) - 


1- a32n-(i) - 


1- a33'?(i) 




-0, 


(B2c) 


^^(1)11 - 


- a4le^l) - 


1- Q;42n(i) - 


1- a43i9(i) 


+ 0:44 "^Ri^i) 


-0, 


(B2d) 




1- Q!5ie(i) - 


1- Q;52n(i) - 


1- Q!53l?(i) 


+ a54^-R(i)|| 


= 0, 


(B2e) 
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Table I: The coefficients aij figuring in tlie equations (B2). 



Qirn , ^ , «^e(o)(l + w) 




e(o)(l + w) 


e(o)(l + w;) 




4iJ 


2H 


on 




n(o) 
AH 


Pe V2 

£((,)(! + w) a2 


Pn V2 

e(o)(l + w) 


i7(2-3/32) 





K "^1^(0) 





-2^6(0) (1 + w) 




-n-(o) 





6e,o,i/(l + w) 


§e(o)(l + w) 


3i?(o) + 3Ke(o)(l + w) 


^i?(o) + 3Ke(o)(l + w) 


^i?(0) + 3Ke(n)(l + w) 



where the coefficients aij{t) are given in Table f. 

fn calculating the coefficients oi, 02 and 03, (198) in the main text, we use that the time derivative of the quotient 
w, defined by (153) is given by 



w 



m{l + w){w- p"^), 



(B3) 



as follows from equation (154c) and the expression (153). Moreover, it is of convenience not to expand the function f5{t) 
defined by (153) since this will considerably complicate the expressions for the coefficients ai, 02 and 03. 

a. Step 1. We first eliminate the quantity from equations (B2). Differentiating equation (B2e) with respect 

to time and eliminating the time derivatives ^(i), fi^^, -d^^ and with the help of equations (B2a)-(B2d), we 

arrive at the equation 



qi) +pi£(i) +P2'^(i) 



Ps-^a) +P4^R(i)\\ = 0, 



where the coefficients pi(t), . . . ,p4(t) are given by 

Pi = 0^51 - 01510^11 - oib2a-2i - asa^Sj - abi<^Ai 

From equation (B4) it follows that 



1 

— i 

PA 



Pi _P2 

£(1) 1^(1) 
Pi Pi 



P3 
Pi 



(B4) 
(B5) 

(B6) 



In this way we have expressed the quantity as a linear combination of the quantities efl^, £(1), n^i-^ and 

Upon replacing ^i?(i)|| given by (B6), in equations (B2), we arrive at the system of equations 



£(1) - 


)- ^lefi) 


+- /3iie(i) 


+ /3i2n-(i) 


+ /3l3^(i) 


= 0, 


(B7a) 




)- '72£fi') 


/32l£(l) 


+ /322n-(i) 


+ /323^9(i) 


= 0, 


(B7b) 




)- qs^fi) 


/33l£(i) 


+ /332»^(i) 


+ /333t?(i) 


-0, 


(B7c) 


'^(1)11 - 


f (74efi') 


■f /34ie(i) 


+ /342^^(l) 


+ /343t?(i) 


= 0, 


(B7d) 


t(i) 


1- Qbefi) 




+ /352n(i) 


+ /353l?(i) 


= 0, 


(B7e) 



where the coefficients qi{t) and Pij{t) are given by 



Pi 



Pii = Oil 



qiPj- 



(B8) 



We now have achieved that the quantity '^R(i)\\ occurs only in equation (B7d). Since we are not interested in the 
non-physical quantity ^R,-^^h, we do not need this equation any more. 
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b. Step 2. We now proceed in the same way as in step 1: eliminating this time the quantity 19(1) from the system 
of equations (B7). Differentiating equation (B7e) with respect to time and eliminating the time derivatives £(1), h(i^ 
and 1^(1) with the help of equations (B7a)-(B7c), we arrive at 



qssfi) + rsfl + sie(i) + S2"-(i) + S3i?(i) = 0, 
where the coefficients r{t) and Si{t) are given by 



From equation (B9) it follows that 



r = 1 + 95 - Pziqi - /352g2 - /^ssgs- 



, _ 95 ..gi r .gi si S2 

C(i) — £(1) £(1) £(1) ^(1)- 

S3 S3 S3 S3 



(B9) 



(BlOa) 
(BlOb) 



(Bll) 



In this way we have expressed the quantity 'd^^ as a linear combination of the quantities ef^, efl^, £(i) and n^^. Upon 
replacing z9(i) given by (Bll) in equations (B7), we arrive at the system of equations 
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r 

S3 



iqi - /?13- 
iqs - /333- 

{qa - /543- 

iqs - ^53- 



t(l) 



- I3l3 
hi ~ /?23 

hi — hs 

Pil — /343 



hi - A 



53" 



£1 
S3 
Si 
S3 
£1 
S3 
Si 
S3 

£1 

■S3 



£(1) 



£(1) + 



/3l2 - ^13 

h2 — hs 

h2 — /333 
/^42 ~ h?, 





= 0, 


(B12a) 




= 0, 


(B12b) 




= 0, 


(B12c) 




= 0, 


(B12d) 




= 0. 


(B12e) 



We have achieved now that the quantities i?(i) and ■^i?(i)|| occur only in equations (B12c) and (B12d), so that these 
equations will not be needed anymore. We are left, in principle, with equations (B12a), (B12b) and (B12e) for the 
three unknown quantities £(i), riji) and £^i), but we first proceed with all five equations. 

c. Step 3. At first sight, the next steps would be to eliminate, successively, the quantities £(i) and rt(i) from 
equation (B12e) with the help of equations (B12a) and (B12b). We then would end up with a fourth-order differential 
equation for the unknown quantity £^') . 

However, it is possible to extract a second-order equation for the gauge-invariant energy density from the equa- 
tions (B12). This will now be shown. Equation (B12e) can be rewritten 



£fi) -I- ai£° ) + a2£fi) = 03 ( n(i) -f 

where the coefficients ai(t), 02 (i) and az{t) are given by 

S3 , r S3 



/351S3 - /353S1 



Oi = 



02 



/352S3 - /353S2^'" 

/352S3 S2 

03 



hi 95 ^5395 ^5395 95 

These are precisely the coefficients (198a)-(198c) of the main text. Furthermore, we find 

his?.-hzsi n,, 



Hence, 



Ha) 

P52S3-h3S2 £(o)(l + U;)' 

hiss - hssi n^o) 



/352S3 - /353S2 



£(0)(1 + w) 



(B13) 



(B14) 



(B15) 



(B16) 



With the help of this expression and (174) we can rewrite equation (B13) in the form (197a). 

The derivation of the expressions (198) from (B14) and the proof of the equality (B15) is straightforward, but 
extremely complicated. We used Maple 13 [45] to perform this algebraic task [46]. 
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3. Evolution Equations for the Contrast Functions 

In this section we derive equations (201). We start off with equation (201b). From (177) and the definitions (200) 
it follows that 

J^)- (B17) 

Differentiating this equation with respect to ct yields 

a,afi, = n<o, - j^) + r.<o, ^ ^ - J^) , (B18) 

where we have used equation (197b). Using equation (154e) and the expression (B17) to eliminate afl-,, we arrive at 
equation (201b) of the main text. 

Finally, we derive equation (201a). Upon substituting the expressions 

efi) = £(0)^6, ^fi) = E{o)Se + £{o}Se, sfi) — ^(a-jS^ + 2i^Q^5^ + e^Qfd^, (-B19) 
into equation (197a), and dividing by Sfo), we find 

Oi = 2 hai, 62 = h ai 1-02, 03 = 03 . (B20) 

£(0) £(0) £(0) £(0) 

where we have also used (B17). These are the coefficients (202) of the main text. 



Appendix C: Gauge- invar iance of the First-order Equations 

If we go over from one synchronous system of reference with coordinates x to another synchronous system of 
reference with coordinates x given by expression (4), we have 

eM;0 + eo;M = 0, (CI) 

as follows from the transformation rule (11) and the conditions (27). From this equation we find, using (28), (30a), 
(30b) and (56) that ^'^{t,x) must be of the form 

where ipix) and x*(a;) are arbitrary functions — of the first-order — of the spatial coordinates x. The fact that the 
gauge function tp does not depend on the time coordinate 2;° = ct anymore, as it did in general coordinates, see (6), 
is a consequence of the choice of synchronous coordinates for the original coordinates as well as for the transformed 
system of reference. 

The energy density perturbation transforms according to (2a), where £(o) is a solution of equation (154c). Similarly, 
the particle number density transforms according to (2b) where n(o) is a solution of equation (154e). Finally, as follows 
from (15), the fluid expansion scalar 6, (20c), transforms as 

e,,,{t,x) ^ 0,,,{t,x) + ^{x)0^„,it), (C3) 

where 6'(o) ~ 3H is a solution of the set (154). 

From (15) with a = p, e, or n and (77) we find for the transformation rule for the first-order perturbations to the 
pressure 

P(i) =Pe£(i) +?«"-(!)■ (C4) 

The transformation rule (16) with the four- velocity u^^ implies 

wfi) - <i) - C^o, (C5) 
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where we have used that ufo^ = S^q, expression (53). From the transformation rule (C5) it follows that u^, transforms 
under transformations (C2) between synchronous coordinates as 

u° ) (t, x) = u° ) {t, x) = 0, (C6a) 

{t, x) - (t, a;) - -L_g^k^x)dki^{x). (C6b) 

We want to determine the transformation rules for 19(1) and Since the quantities ^R, (42), and (47), are 

both non-scalars under general space-time transformations, the transformation rule (15) is not applicable to determine 
the transformation of their first-order perturbations under infinitesimal space-time transformations f , (4). 

Since u^^jy satisfies equation (141e), and since M^j^jy transforms according to (C6), and since we know that u'^^y 
satisfies equation (141e) with hats, one may verify, using (C4), that 

{,^,^{t,x)^d,,,{t,x)~^^^, (C7) 

satisfies equation (161d) with hatted quantities. The quantity i9(i) is defined in analogy to -d^^ in (74) 

^(1, - («a,||)|fe- (C8) 

Apparently, i9(i) transforms according to (C7) under arbitrary infinitesimal space-time transformations between syn- 
chronous coordinates. Similarly, one may verify that 

%,,\\{t,x)^^R,,,\\{t,x) + m{t) |^^^-i3i?,o,(i)^(a;)j , (C9) 

satisfies equation (161a). Apparently, expression (C9) is the transformation rule for under arbitrary infinitesimal 

space-time transformations between synchronous coordinates. An alternative way to find the results (C7) and (C9) 
is to write d^i) = — f and ^R(i)\\ = "^^(1)11 ~5i where / and g are unknown functions, to substitute, thereupon, ijji) 
and into equations (161d) and (161a), and to determine / and g such that the old equations (161d) and (161a) 

reappear. In fact, our method to define i9(i), (C7), and (C9), is nothing but a shortcut to this procedure. 

It may now easily be verified by substitution that if £(1), n(i), ^(i), ?9(i), and are solutions of the systems 

(152) and (195), then the quantities £(1,, (2a), n(i), (2b), 6^1-,, (C3), (C7), and '^R(i)\\, (C9), are, for an arbitrary 
function ip^x), also solutions of these systems. In other words, the systems (152) and (195) are gauge-invariant 
under gauge transformations between synchronous coordinates. The solutions £(i), n^i^, and however, 

contain an arbitrary function tpix) and are, therefore, gauge dependent. 



Appendix D: Horizon Size after Decoupling of Matter and Radiation 

As is well-known, a density perturbation can only grow if all its particles are in causal contact with each other, so 
that gravity can act in such a way that a density perturbation may eventually collapse. In this appendix we calculate 
the horizon size at time tmat between the decoupling time tdec and the present time tp. The horizon size at time t is 
given by 

r* dt' 

dn{t)^ca{t) — (Dl) 
Jo a{t') 

Using (269), we get 

du{t) = 3ct. (D2) 

For the horizon size at time tmati we find 

NO 3ctp 1.260 X 10^ , 

aH(tmat) = 3ct,nat = ^ = ^ kpC, (D3) 

where we have used (231), (232) and (269). At decoupling, 2;(tdoc) = 1091, the horizon size is c?H(^dcc) — 349 kpc. 
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Appendix E: Derivation of the Relativistic Standard Equation for Density Perturbations 

In this appendix equations (355) and (362) of the main text are derived, for a flat flrw universe, ^R^o) — 0, with 
vanishing cosmological constant, A = 0, using the background equations (154) and the hnearized Einstein equations 
and conservation laws for scalar perturbations (195). 

From (B3) it follows that w is constant if and only if w = for all times. Using (182) it is found for constant 
w that p„ = and = w, i.e. the pressure does not depend on the particle number density. Consequently, in the 
derivation of equations (355) and (362) the equations (154e) for n^g^{t) and (195b) for n^i^{t,x) are not needed. In 
this case, the equation of state is given by 

p — we. (El) 
To derive the standard equations (355) and (362), it is required that Wj^jy ~ 0, implying that 

d,,,{t,x)^0, ^(a;) = V, (E2) 

where we have used (218). The first of (E2) implies, using, equation (195c) that 

V^Pa) = 0, (E3) 

i.e. pressure gradients should vanish in order to derive the relativistic standard equation (E6). Substituting e^if — £(o)(5 
into equation (195a) and eliminating e^Q) with the help of equation (154c), it is found that 

'^+^Ko,'5+ri?a,||) -0. (E4) 

Differentiating equation (E4) with respect to — ct and using equations (154) and (195d), yields 

^ + 1(1 + w)H5 - 1(1 + w)^Ke,„^S - 1(1 + w){l - 3w) ^i?,,,,, = 0. (E5) 

Eliminating from equation (E5) with the help of equation (E4), yields the standard equation for large-scale 

perturbations in a flat FLRW universe: 

6 + 2Hd-^Ke^o)S{l + w){l + 3w)^0. (E6) 

This equation has been derived by Weinberg [35], equation (15.10.57) and Peebles [9], equation (86.11). For w — ^ 
(the radiation-dominated era) equation (355) is found, whereas for w — (the era after decoupling of matter and 
radiation) equation (362) applies. 

Using that the general solution of the background equations (154) for '^R(o) = 0, A = and constant w is given by 



m - ^(T^(^^)"' - ^(^o) i^^ ] , (E7a) 



e(o,(i) = ^ .^ ict)-^ = eM (I) , (E7b) 

SkII + w)^ \to J 

we find for the general solution of equation (E6), with t = t/to, 

S{t) = Sir(2+6»)/(3+3-) _ 3(1 + w)H{to)ijT-\ (E8) 

where Ei is an arbitrary integration constant. 

The only surviving gauge mode is, using (2a), (154c) and 5 = En^/e^o) [cf. (200)], 



e,r„{t) = Ve(o,(i), W(0 = ^'^'"^^^ = '/^^ = -3(1 -t- w)H{t)^. (E9) 



It follows from (E7a) and (E9) that the second term in the right-hand side of (E8) is a gauge mode and tp is the gauge 
constant. For w = ^ we find (354) with gauge mode (356) and for w = we get (361) with gauge mode (363). The 
solution (E8) is exactly equal to the result found by Peebles ]9], §86, expression (86.12). 

Finally, we note that in the derivation of (355) [i.e. (E6) with w = |] Peebles uses = (in his notation: 6 = 0), 
see (E2). In this case, Peebles' method yields a physical mode S (x t and a gauge mode S oc t~^. For ^ Peebles 
finds a physical mode 8 cx t^I"^ . However, in the treatise presented in this article both physical modes (5 cx r and 
(5 cx r^/^ ]see (258)] follow from one second-order differential equation (245a). Moreover, in our treatise neither the 
pressure gradient Vp^',, nor the velocity divergence = (Mj^jj||)|fe are neglected in our evolution equations (201). 
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Appendix F: Symbols and their Meaning 

Table II: Symbols and their meaning of all quantities, except for those 
occurring in the appendices. 



(0) 

(1) 

(2) 

;A 

Ifc 



Symbol 
V/ 

iyfY 

V • u 

(V A u)i 

vV 

do 
ft 

hat: " 
dot: ■ 
prime: ' 
tilde: ~ 

super-index: ^' 
super-index: 
sub-index: 
sub-index: 
sub-index: 
sub-index: 
sub-index: 
sub-index: 
sub-index: y 
sub-index: ± 
sub-index: , 

/3 

7 
5, 

Sn 
Sp 

St 

£ 

e 

e 

K 

A 
A 
M 

Mr 
Mm 

TT 
W 

P(l) 
O" 

.Q...;3 
a 

as 

ai, 02, as, o-i 
bi,b2, bs 



Meaning 

{dif,d2f,d3f) 



Ct 



Ct 



V-(V/)=rVN. 
derivative with respect to x" 
derivative with respect to 
computed with respect to x 
derivative with respect to x" 
derivative with respect to r 
computed with respect to three-metric gij 
gauge-invariant 

contravariant derivative with respect to x*^: (^'* = g'^^^.^(\j 

background quantity 

perturbation of first-order 

perturbation of second-order 

covariant derivative with respect to 

covariant derivative with respect to x'^ 

derivative with respect to x'^ 

longitudinal part of a vector or tensor 

perpendicular part of vector or tensor 

transverse and traceless part of a tensor 

c^^ times speed of sound: y^P(o) /£(o) 

connection coefficients 

arbitrary function 

energy contrast function 

Kronecker delta 

particle number contrast function 

pressure contrast function 

(matter) temperature contrast function 

background radiation temperature contrast function 

energy density 

energy per particle 

Levi-Civita tensor, ei'^"^ = -1-1 

potentials due to relativistic density perturbations 

bookkeeping parameter equal to 1 

expansion scalar in four-space 

expansion scalar in three-space 

SttG/c"* 

time derivative metric coefficients 
cosmological constant 
wavelength, physical scale at time tp 
thermodynamic potential 
reduced wave-number (radiation) 
reduced wave-number (matter) 
first-order space-time translation 
arbitrary function 
FLRW-coordinate 

arbitrary scalar 
abbreviation for n^J^ 
dimensionless time 

potentials due to relativistic density perturbations 
Newtonian potential 
arbitrary three- vector 
first-order time translation 
components in the Friedmann equation 
arbitrary scalar 
arbitrary tensor 
arbitrary rank two tensor 
scale factor or radius of universe 
black body constant 
coefficients in perturbation equations 
coefficients in perturbation equations 



■n(Q)£nJ[e(o)il + w)] 



Reference Equation 

(125)-(128) 
(125), (126) 
(128) 
(151) 



(2) 
(28) 

(253), (298) 
(51b), (62) 

(3) 

(111a) 

22) 

22) 
(50) 
(7) 
(33) 

(153) 

(29) 

(113a) 

(200) 

(53) 

(200) 

(206) 

(205), (286) 
(205), (336) 
(20a) 
(164) 
128) 

(111a), (207) 

(50) 

(20c) 

(47) 

(38) 

(28) 

(36) 

(248) 

(167) 

(254) 

(295) 

(113a) 
(52) 
(229) 
(10a) 

(176), (177) 
(251), (296) 
(111a), (207) 
(225), (224) 
(C2) 

(2), (4), (6) 

(155) 

(18) 

(10c) 
51b) 
(236) 
(198) 

(202) 
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X 



Table II: (continued) 



~c speed ot light 

ds'^ line element in four-space 

E energy within volume V 

G Newton's gravitation constant 

gi_ii, metric tensor 

Qi^i, metric with respect to x'^ 

Qij time-independent metric of three-space 

H c-^n 

H Hubble function: H = (da/dt)/a 

hij minus first-order perturbation of metric 

k k = —1 (open), (flat), +1 (closed) 

ks Boltzmann's constant 

£^ Lie derivative with respect to 

m mass of particle of cosmological fluid 

M© solar mass, 1.98892 x 10^° kg 

niH baryonic (proton) or cdm mass 

A'^ number of particles within volume V 

particle density four-flow 

n particle number density 

p pressure 

Pe,Pn partial derivatives of pressure 
PnnjPenyPee partial derivatives of pressure 

q circular wave number: 2tt/X 

r FLRW-coordinate 

Ricci scalar in three-space 

Rfi^ Ricci tensor in four-space 

^Rij Ricci tensor in three-space 

S entropy within a volume V 

s entropy per particle s = S/N 

s^j gauge-invariant entropy perturbation 

T absolute temperature 

T-y photon temperature 

T^" energy momentum tensor 

Tj^j'j gauge-invariant temperature perturbation 

t cosmological time 

to initial cosmological time 

tdcc decoupling time 

tcq matter-radiation equality (onset of matter-dominated era 

tmat initial time between decoupling and the present 

tp present cosmological time (13.7 Gyr) 

trad onset of radiation-dominated era 

If^ cosmological four- velocity (^C/^ = (20) 

c-'C/" — , 

U spatial velocity (212) 

u c~^U (124 
arbitrary four- vector (lOb 

Vs speed of sound (273) 

w pressure divided by energy density (153) 

X space-time point x'^ = [ct, x) — 

x^ locally transformed coordinates (4) 
spatial point x — {x^,x'^,x^) 



z(t) redshift at time t (232) 
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